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Die numerische Auflésung des zweiten und des dritten Randwertproblems 
mit Hilfe der Differenzenmethode ist gegeniiber dem ersten Randwertproblem > 
bis jetzt stiefmiitterlich behandelt worden. Wohl hat L. F. RicHARDSON [16] 
Beispiele dargeboten. Auf ein numerisches Verfahren sind aber erst F. WOLF [24] 
im Falle eines Stufenrandes und R. V. SOUTHWELL, G. VAISEY [20], [21] bei 
gekriimmtem Rand eingegangen. Die letztere Methode kniipft stark an mecha- 
nische Vorstellungen an, die jedoch beim gemischten Randwertproblem kaum 
mehr verwendet werden kénnen. Neuerdings hat sich L. Fox [7], [8] mit den 
Differenzengleichungen, die sich in den Gitterrandpunkten ergeben, besonders 
beschaftigt. Doch blieb das Problem der Approximation bei gekriimmtem 
Rande noch offen. : 

In § 1 studieren wir zunachst den Ubergang zu Differenzengleichungen von 
neuem, wobei wir der Approximation am Rande besondere Aufmerksamkeit: 
schenken. Allerdings sehen wir uns gezwungen, Ecken und Spitzen am Rande 
auszuschliessen, sofern nicht gerade der einfache Fall eines Stufenrandes vor- 
hegt. Die eigenartigen Schwierigkeiten, welche eine Unstetigkeit der Tangente 
erzeugt, kommen in Abschnitt 8 zur Sprache. ; 

Wir geben auch in § 1 hinreichende Bedingungen fiir die Existenz und Ein- 
deutigkeit von Gitterlésungen, das heisst der Lésungen des dem gegebenen 
Problem zugeordneten Systems von Differenzengleichungen. Wir erwahnen 
hier nur nebenbei, dass man, von unseren Ergebnissen ausgehend, sehr wahr- 
scheinlich entsprechende Existenz- und Unitatssatze fiir Differentialgleichun- - 
gen beweisen kénnte, wie solche von R. Courant, K. FriepRicHs, H. Lewy [5] 
im Falle von Randbedingungen erster Art hergeleitet worden sind. 

Uber den beim Differenzenverfahren begangenen A bbrechfehler (truncation: 
error), das heisst iiber die Differenz zwischen der Gitterlésung und der gesuch- 
ten kontinuierlichen Lésung, hat S. GERSCHGORIN[10] im Falle des ersten 
Randwertproblems gezeigt, dass dieser Fehler bei Verfeinerung des Gitters’ 
mindestens proportional zur Gitterseite gegen Null strebt. L. CoLLATz [1] hat 
anschliessend dargetan, dass eine bessere Anpassung der Randwerte durch 
Interpolation die Konvergenz wesentlich verbessert, so namlich, dass der Fehler 
mit h® gegen Null strebt. Ein solches Ergebnis war schon friiher erwartet wor- 
den, wie Andeutungen von L. F. RICHARDSON [16] und von C. RUNGE [17] 
zeigen, 

In § 2 werden wir eine Fehlerabschatzung fiir den Fall von Randbedingun- 
gen dritter Art sowie auch fiir gemischte Randbedingungen vornehmen. Wir 
modifizieren und erweitern hiezu die von S. GERSCHGORIN [10] ersonnene Majo- 
rantenmethode, Der Einfachheit halber beschriinken wir uns auf Differential- 
gleichungen zweiter Ordnung mit zwei unabhangigen Variablen, wobei wir 
Eigenwertprobleme ausschliessen. In Satz 1 (§ 2, Abschnitt 14) sagen wir aus, 
dass bei dem von uns vorgeschlagenen numerischen Verfahren der Abbrech- 
fehler proportional zur Gitterseite h gegen Null strebt. Anschliessend daran 


- 
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diskutieren wir die Massnahmen, die man zur Konvergenzbeschleunigung zu 
treffen hat. Unser Vorgehen ist an allerlei einschrankende Voraussetzungen 
gebunden, welche in erster Linie die Existenz einer eindeutigen Losung mit 
gewissen Stetigkeitseigenschaften sicherzustellen haben. Immerhin ist zu er- 
warten, dass das Differenzenverfahren auch in allgemeineren Fallen brauchbare 
“numerische Resultate liefert, wie praktische Erfahrungen zeigen. Doch wird 
der Konvergenzgrad im allgemeinen geringer sein. 
Schliesslich gehen wir in § 3 auf die numerische Auflésung des zageordneten 
Systems von Differenzengleichungen genauer ein. Wir studieren ein alternie- 
rendes Verfahren, das als eine Erweiterung des Southwellschen Relaxationsver- 
fahrens aufgefasst werden kann. In Satz 2 (Abschnitt 18) zeigen wir, dass dieses 
Approximationsverfahren fiir jedes feste Gitter konvergiert. In Abschnitt 19 
-folgt eine Diskussion der Konvergenzgeschwindigkeit. Beim Beweis von Satz 2 
machen wir ausgiebig von den Greenschen Zahlen Gebrauch, welche wir weiter 
-vorne, in Abschnitt 6, einfiihren. Die Arbeit schliesst mit dem Studium der 
unyermeidlichen Abrundungsfehler. 


§ 1. Das der Differentialgleichung und den Randbedingungen 
entsprechende System von Differenzengleichungen 


1. Es bezeichne G ein beschranktes, einfach oder mehrfach zusammenhan- 
~ gendes Gebiet der reellen (x, y)-Ebene. Der Rand C des Gebietes soll frei von 
Doppelpunkten, Ecken und Spitzen sein, also iiberall eine stetige Tangente 
_besitzen. ; 
Gesucht ist eine Funktion v = v(x, y), welche in G der normierten!) ellipti- 
- schen Differentialgleichung 


0?u 0?u Ov Ov 
Ne iat es? she aes z= 1 
geniigt und auf C die Randbedingung 
Ov 
— - —— — ie 
n(s) 22 — hls) 0 = #(s) (1. 2) 


bei nach innen gerichteter Normale erfillt. 

Wir setzen voraus, dass die Koeffizienten a, 0, c, gin (1.1) solche beschrankte 
Funktionen von x und y und A(s), A(s), /(s) in (1. 2) solche beschrankte Funk- 
tionen eines auf C erklarten Parameters sind, dass das Problem eine eindeutig 
bestimmte Lésung v(x, y) besitzt, die im Innern und auf dem Rand von G samt 


1) Uber die Transformation der allgemeinen elliptischen Differentialgleichung auf die Normal- 
- form siehe etwa CouRANT-HILBERT, Methoden der mathematischen Physik, Bd. 2 (Springer, Berlin 


1937), S. 128 ff. 
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ihren partiellen Ableitungen bis zur vierten Ordnung stetig ist+). Die verlangte 
Stetigkeitseigenschaften sollen nicht verlorengehen, wenn g und f(s) durch ein 
gemeinsame, von Null verschiedene Konstante ersetzt werden?). Unter allen 
Umstanden soll im folgenden stets 


e=c(z, vy) S0 (1. | 


gelten. Ferner miissen wir fiir die Randbedingung die folgenden Annahmen 
machen: 

Der Rand C zerfalle in zwei Teile C, und C,, wobei C, und C, aus endlich 
vielen Teilbégen bestehen kénnen. Dann sei 


k(s} =leaut Cy, “R(s) = Oran G, (1. 4) 
und 
h(s) 20 aufC,, h(s)=1 aufC,. (1. 5) 


Se ram 


Sollte die Menge der Randpunkte C, leer und zugleich in allen inneren 
Punkten c(x, y) = 0 sein, so verlangen wir iiber (1. 5) hinaus, dass auf wenig- 
stens einem Teilbogen von C 

h(s) > 0 (1. 6) 
gilt. 

In dem soeben formulierten Problem sind die folgenden Spezialfalle mit- 
enthalten: 

a) Das Randwertproblem erster Art, wenn namlich C, eine Nullmenge ist. 

b) Das Randwertproblem dritter Art, wenn C, eine Nullmenge ist. 

c) Das aus Randbedingungen erster und zweiter Art hervorgehende ge- 

mischte Problem. Hier ist langs C, die Normalableitung 0v/0n und langs 
C,, das keine Nullmenge ist, der Wert von v vorgegeben. 

Das zweite Randwertproblem ist nur dann mitenthalten, wenn c(x, y) nicht 
in allen Punkten von G verschwindet. Denn sonst muss fiir die Existenz einer 
Lésung bekanntlich 


ati Aa a aaa nil lla aan ale tag eA 


vorausgesetzt werden. 


1) Unter welchen Bedingungen fiir die Koeffizienten eine eindeutige Lésung mit den verlangten — 
Stetigkeitseigenschaften existiert, sei hier nicht untersucht. Die Literatur iiber Existenzsatze fiir — 
das dritte Randwertproblem findet man zusammengestellt im Artikel von L. LicHTENSTEIN 
Enzyklopddie der mathematischen Wissenschaften, Bd. 2, 8. 2, S. 1803-1305, oder im Artikel von 
W. STERNBERG in Pascars Repertorium der hiheren Mathematik, Bd. 1, 3, 2. Aufl. (Teubner, Leipzig 
und Berlin 1929), S, 1189. Es ist dort iiblich, die auch von uns gemachten Voraussetzungen (1. 3) 
und (1. 6) zugrunde zu legen. Doch wird die Stetigkeit der héheren als zweiten Ableitungen nicht 
verlangt. Auf der anderen Seite ist unser Problem mit Riicksicht auf die Anwendungen allgemeiner 
formuliert, indem k(s) auch den Wert Null annehmen darf. 7 

*) Da die linken Seiten in (1. 1) und (1. 2) homogen in v sind, so ist die Wahl der Konstanten 
unwesentlich. In (13. 1) wird sie gleich —1 gesetzt. 
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2. Um (1. 1) mit der Randbedingung (1. 2) numerisch auflésen zu k6nnen, 
wenden wir das Differenzenverfahren an. Zu diesem Zweck legen wir in die 
(x, y)-Ebene ein quadratisches Punktgitter parallel zu den Koordinatenachsen!), 
Die Seitenlange, das heisst der Abstand zweier benachbarter Gitterpunkte, sei 
mit / bezeichnet. Den Rand C von G approximieren wir auf die folgende Weise 
durch einen von den Gitterpunkten gebildeten Streckenzug C’ (Figur 1): 


Fig. 1 
Approximation des Gebietes G durch das Gittergebiet G’. 


Von jeder Gitterseite, die von C geschnitten wird, ordne man den dem 
Schnittpunkt am nachsten liegenden Gitterpunkt dem Rande C’ zu. Im Zwei- 
felsfall bleibt die Wahl frei. Von einem so bestimmten Gitterrandpunkt kann 
man auf C ein Lot fallen, welches < h/2 ist. Dies folgt unmittelbar aus der Vor- 
aussetzung, dass C iiberall eine stetige Tangente besitze. # sei dabei so klein 
_gewahlt, dass zu jedem Punkt auf C’ nur ein einziges solches Lot existiert. Auf 
diese Weise wird jedem Punkt Q von C’ eindeutig ein Punkt Q von C, namlich 
der Fusspunkt des Lotes, zugeordnet. Mit Cj und Cj bezeichnen wir die Mengen 
derjenigen Punkte von C’, denen ein Randpunkt von C, bzw. C, zugeordnet ist. 
Fiir hinreichend kleines / sind Ci und C; nicht leer, wenn das Entsprechende 
fiir C, bzw. C, gilt. 

Mit G’ bezeichnen wir das von C’ eingeschlossene Gittergebiet. Jeder innere 
Gitterpunkt ist dadurch ausgezeichnet, dass er die maximale Zahl von vier 
Nachbarpunkten besitzt, die alle zum «abgeschlossenen» Gittergebiet G’ ge- 
héren. Wir diirfen erwarten, dass / bereits so klein ist, dass die innern Gitter- 
punkte von G’ eine «zusammenhangende» Punktmenge bilden, und dass auch 
keine «Engpiisse» existieren. In G’ bilden wir nun Differenzenquotienten, welche 


1) Gelegentlich eignen sich andere als quadratische Gitter besser. Wir verweisen hiezu aut 
R. V. SouTHWELL [21]. 
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die in (1. 1) und (1. 2) vorkommenden Differentialquotienten approximieren. 
Eine Gitterfunktion u(x, y) nehme in einem beliebigen Gitterpunkt /% den 
Wert w#, und in den vier Nachbarpunkten P, P,, P;, P, die Werte um, ..., M4 an 
(Figur 1). Einige dieser Nachbarpunkte kénnen auch Randpunkte sein. Unter 
den partiellen Differenzenquotienten erster und zweiter Ordnung verstehen wir 
dann wie iiblich die Ausdriicke ; 


a5 aE oie ~ “sin der x-Richtung, | 
7 (2. 1) 

Pas 3G Fea! : ; 
= os a a + “4 in der y-Richtung. ; 


Ersetzt man in (1. 1) die Differentialquotienten durch die soeben eingefiihrten. 
Differenzenquotienten, so lauten die zu (1. 1) analogen Differenzengleichungen 


Uy — Us 


2h = +2 


Ug — 


2h 


1 
sr (Ur + ue + Mg + Uy — 4 Ug) + a “t+ Cum=e, 


j 
wobei die Werte von a, b, c, gim Punkte P, zu bestimmen sind. Daraus gewin- 
nen wir durch Umformung : 

: 


ah bh ah’ 
(1 + =) Uy + (1 + =) Us + (1 = 1 Us | ; 
Me (2. 2) 
+ (1 - =) u,— (4—ch*)u,=h'¢, | ! 
oder abkiirzend die Differenzengleichungen 
Oy Uy + My Ug + Hg Ug + Oy Uy — BUy=y. (2533 . 


Hierin sind die Koeffizienten a, %, %3, «4, 8, y gewisse Gitterfunktionen, welche 
wegen (2. 2) und (1. 3) den folgenden Bedingungen geniigen: 


p=4. (2. 5) 
Es sei ausserdem / so klein angenommen, dass iiberall in G’ 


ag 0 (¢ = 1, 2,3, 4), Gay 
gilt. 

3. Um die Normalableitung passend approximieren zu kénnen, sei Q ein 
beliebiger Gitterpunkt von C’ und n = n(s) die durch Q gehende Kurvennor- 
male von C (Figur 1). Durch sie wird, wie in Abschnitt 2 dargelegt worden ist, 
dem Gitterpunkt Q ein Punkt Q der Kurve C zugeordnet. Folgt man der Kur- 
vennormalen 7 ins Innere von G’, so schneidet sie die nachstfolgende Gitter- 
seite in einem gewissen Punkt S. S ist im allgemeinen kein Gitterpunkt. Die 
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Endpunkte der Gitterseite, auf der S hegt, seien mit P’ und P” bezeichnet. Bei 
hinreichend kleiner Gitterseite # ist mindestens einer dieser beiden Gitter- 
punkte ein innerer Punkt von G’. 

Weiterhin sei der Abstand des Punktes S von Q mit o h bezeichnet, wobei 
fur o offenbar 


isoss2 (3. 1) 


gilt. Aus den Werten uw’ und w”, welche u(x, y) in den Punkten P’ und P’ 
annimmt, werde durch lineare Interpolation der Wert u in S gebildet, so dass 
also uw die Gestalt 


t=) Wag", AptAg=1, 14,20, 4,20 (3. 2) 


hat (vgl. Figur 3 in Abschnitt 10). Mit U bezeichnen wir den Wert von u(x, y) 
in Q. Dann werde 0v/dn im Punkte Q durch den Differenzenquotienten 


u— U 


ar (3.3) 
im Punkte Q approximiert. 
Die Randbedingung (1. 2) wird daher im Punkte Q durch 
u—U 
k(s) ei ska h(s) U = f(s) (3. 4) 
angenadhert, woftir wir nach Multiplikation mit o # auch 
k(s) wu —[k(s) + ch h(s)| U = oh f(s) (3. 5) 
oder kurz 
Ku-HU=F (3. 6) 


schreiben kénnen. Darin bezeichnen K, H, F gewisse Gitterfunktionen auf dem 
Rand C’. Die zu ihrer Bestimmung im Randpunkt Q notwendigen Werte von 
R(s), A(s), f(s) sind in dem zugeordneten Punkt Q zu berechnen. 
An Stelle von (1. 4) tritt die Annahme 


Kz laufCl, K=0 aul (3. 7) 


und an Stelle von (1. 5) 

Hew oh is) 2 autC,, H=oh > 0 aut Co, (3. 8) 
wobei fiir den Fall, dass Ci eine Nullmenge ist, / so klein vorauszusetzen ist, 
dass fiir wenigstens einen Punkt von C’ ; 


Jel > il (3. 9) 
gilt. 
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4. Die Gleichungen (2. 3) und (3. 6) fassen wir zu einem System zusammen, 
wobei wir gleichzeitig fiir die Linearformen linker Hand abkiirzende Bezeich- 
nungen einfithren: 


L(U)=Ku—HU=F. 


(4. 1) 


oe 


Hierin verwenden wir grosse Buchstaben nur fiir solche Gréssen, die sich 
auf den Rand C’, kleine Buchstaben aber fiir Gréssen, die sich Morgue auf 
das Innere von G! beziehen. 

Bezeichnet N die Anzahl aller Gitespanite von G’ mit incenhiey des 
Randes C’, so bilden die Gleichungen (4. 1) ein System von genau WN linearen 
Gleichungen mit den N Unbekannten u(x, y). 

Dieses System besitzt unter den Bedingungen (2. 4), (2.5), (2. 6), (3.7), 
(3. 8), (3. 9) eine Reihe wichtiger Eigenschaften, die wir in den folgenden Hilfs- 
satzen zusammentfassen : 

Hilfssatz 1. Machen wir im Gleichungssystem (4. 1) zusdtzlich die Annahme 


y20, F20, (4. 2) 


so gilt fiir eine allfdllig existierende Lisung von (4.1) im Innern und auf dem 
Rand von G’ 


w <0. (4. 3) 


Beweis. Angenommen, es existiere eine Lésungsfunktion u(x, y) mit wenig- 
stens einem positiven Wert. Dann bezeichne u*-> 0 den maximalen Wert von 
u in G’, Betrachten wir zuerst die Konsequenzen fiir den Fall, dass u* einem 
Punkt im Innern von G’ zugeordnet ware. Fiir diesen Punkt miisste dann nach 
(4. 2) 

Ay Uy + Ky Ug + Hg Ug + %y Uy — Bu* =O . 
: 
gelten, wobei sich unter den Nachbarpunkten auch Randpunkte befinden kén- : 
nen. Daraus ergibe sich fiir einen beliebigen der vier Nachbarpunkte, etwa fiir 
u, wegen (2. 4), (2. 5), (2. 6) 


Oy Uy = (B — Oy — Oy — Oo) U* = (B— 4+ «,) ut 
oder 


4 
Hy 


m 2 (1+ 2 Jue Sur, (4. 4) 
Sollte im betrachteten Punkt 6 > 4 sein, so kénnte (4. 4) nur mit dem Grésser- 
zeichen bestehen. Damit hatten wir aber wegen der Maximumeigenschaft von 
u* bereits einen Widerspruch. Fiir 6 = 4 wiirde man uw, = vu. = We u = ue 
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erhalten. Man kénnte dann den gleichen Schluss, der zu (4. 4) gefiihrt hat, mit 
den im Innern von G’ liegenden Nachbarpunkten fortsetzen. Da G’ zusammen- 
hangend ist, so werden davon alle Gitterpunkte erfasst. Einem Widerspruch 
k6nnten wir vorerst nur in dem Sonderfall entgehen, wo f iiberall gleich 4 ist. 
_Hierbei hatte «(x, y) im Innern und am Rande den konstanten Wert u*. _ 
Aus dem Gesagten diirfen wir schliessen, dass der maximale Wert u* in 
wenigstens einem Punkt des Randes angenommen wiirde. Ein solcher Punkt 
kann nicht zu Cj gehoren. Denn dort miisste, wegen K = 0, H > 0, F => 0 und 
(4. 1), U <0 sein. Es wiirde sich also um einen Punkt von Cj handeln. Daher 
ware = 1, und man kénnte (4. 1) die Ungleichung 


u 2H u* (4. 5) 


-entnehmen, wo u wegen (3. 2) im allgemeinen einen Mittelwert zweier Funk- 
tionswerte wu’ und uw” bezeichnet. Nun ist wegen (3. 8) H = 1, weshalb (4. 5) 
in allen Fallen zum Widerspruch fiihrt, wo nicht gerade H=1 und 
u=u' =u" = u* ist. In diesem Sonderfall aber gabe es einen inneren Punkt 
mit dem maximalen Funktionswert u*, so dass nach den aus (4. 4) gezogenen 
Konsequenzen tiberall 6B = 4 und u(x, y) = u* gelten miisste. Daher ware C3 
eine Nullmenge. Nach (3.9) gabe es dann wenigstens einen Randpunkt mit 
H > 1. Doch ware auch fiir einen solchen Punkt der LOS ae in der Un- 
gleichung (4. 5) offenkundig. 

5. Hilfssatz 2. Verschwinden alle Werte von y in G' und alle Werte von F 
auf C’, ist also das Gleichungssystem (4. 1) homogen, so besitzt es nur die triviale 

Losung u(x, y) = 0. 

Beweis. u(x, y) =0 ist sicher eine Lésung des Gecroecnen Gleichungs- 
systems. Gabe es noch eine von ihr verschiedene Lésung u*(x, y), so miusste 
wegen Hilfssatz 1 w*(x, y) <0 sein. Da aber mit u*(x, y) auch — u*(x, y) eine 
Lésung ware, so miisste u*(x, y) 20 zugleich mit u*(x, y) [0 gelten, was 
u* = 0 zur Folge hatte. 

Hilfssatz 3. Das Gleichungssystem (4.1) besitzt stets eine eindeutig be- 
stimmte Losung. 

Beweis. Aus dem Hilfssatz 2 folgt, dass die Determinante des Systems (4. 1) 
von Null verschieden sein muss. Dies ist aber eine hinreichende Bedingung fiir 
die Existenz einer eindeutig bestimmten Lésung. 

Hilfssatz 4. Es seien u(x, y) und w(x, y) zwei Gitterfunktionen in G’, fiir 
welche 

Iw) Siu), LW) SLO) (65.1) 


vorausgesetzt ist. Dann gilt im Innern und auf dem Rand von G' 


W2U. (52) 
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Beweis. Man fiihre die Gitterfunktion z= «—w ein. Aus (5.1) folgt — 
I(z) = 0, L(Z) = 0, so dass z die Bedingungen (4. 2) des Hilfssatzes 1 erfiillt. . 
Somit ist z < 0, w = u, q.e. d. 

Hilfssatz 5. Es seien u(x, y) und w(x, y) zwei Gitterfunktionen in G’, welche — 
die Bedingungen 


Iw) <—|Hu)|, LW) S-|LW)| (5.3) 
erfiillen. Dann gilt im Innern und auf dem Rande von G' 
w 2 |u|. - .(6.4)m 


Beweis. Dem Hilfssatz 4 entnimmt man, dass bei Verkleinerung von /(w) und — 
von L(U) die Gitterfunktion u(x, y) wachst. Nun folgt aus (5. 3) /(w) S/(w), 
L(W) <L(U) und damit w =w. Anderseits ergibt sich aus (5. 3) durch © 
Zeichenumkehr /(—w) = |/(u)|, L(—W) =|L(U)| und daher /(—w) 2 J(u) , 
L(—W) = L(U),, was mit Hilfssatz 4 zu —w<u oder zu w= —u fihrt. 7 
Die Ergebnisse w = u und w = —w sind in der Ungleichung (5. 4) zusammen- 
gefasst. 

6. Fiir spatere Konvergenzbetrachtungen fiihren wir im Gittergebiet G’ — 
spezielle Gitterfunktionen g,(x%, y) ein. G’ besitze » Randpunkte, die wir in 


irgendeiner Reihenfolge mit Q, (v = 1, 2,..., 2) bezeichnen. Dann werde g,(x, y) 
folgendermassen definiert : 
a) Es sei y eine der Zahlen 1, 2,...,. Dann nehme g,(x, y) in Q, den Wert 1, 


in allen tibrigen Randpunkten den Wert Null an. 

b) g,(x, y) geniige in jedem innern Punkt von G’ der Differenzengleichung 
l(u) = 0 mit den in (2. 4), (2.5) und (2. 6) gemachten Annahmen iiber 
die Koeffizienten. 

Wegen a geniigt g,(x, y) einer (inhomogenen) Randbedingung erster Art. 
Daher ist g,(x, y) nach Hilfssatz 3 eine eindeutig bestimmte Gitterfunktion. 
Die obige Definition erinnert an diejenige einer Greenschen Funktion. Wir 
wollen deshalb die Werte von g,(x, v), (vy = 1, 2,..., 2), die Greenschen Zahlen 
des Gittergebietes G’ nennen?),. 

Fiir alle Punkte im Innern von G’ leiten wir die Ungleichung 


0< g(x,y) <1 (6. 1) 
her. Zu diesem Zweck greifen wir einen beliebigen innern Punkt P, heraus. Dann 


besagt die Gleichung /(v) = 0, der g,(x, y) = g,(P) definitionsgemass geniigt, 
nach einfacher Umformung: 


1) Fiir den Spezialfall der Potentialgleichung wurden «Greensche Zahlen» schon mehrfach ein- 
gefulhrt. H, Lirpmann [11] und D. Mosxovirz [13] geben ihre Werte fiir einfache quadratische und 
rechteckige Bereiche an. F. Woxr [24] spricht von «Gewichten» und verwendet sie fiir seinen Kon- 
vergenzbeweis des Liebmannschen Verfahrens. 
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g,(B) = ; , el &(Pi) + % &y(P2) + Os Sy(Py) + oy (Pa) ; (6. 2) 
Da die vier Koeffizienten « nach (2.6) und (2. 4) als positive Zahlen mit der 
Summe 4 vorausgesetzt sind, so bedeutet der zweite Faktor rechts in (6. 2) ein 
arithmetisches Mittel mit den Gewichten «. Diese Grésse liegt daher im Innern 
des von den Zahlen g,(FR), ..., g,(P,) gebildeten Intervalls. Da ferner g,(x, y) 
nach Definition am Rande nur der Werte 0 und 1 fahig ist, so liegen siimtliche 
Mittelwerte im Innern des Intervalls (0, 1). An dieser Eigenschaft andert der 
andere Faktor 4/8 wegen (2. 5) nichts, womit (6. 1) bewiesen ist. 

Wir betrachten jetzt alle Gitterfunktionen g,(x, y) (v= 1, 2,..., ”). Ihre 
Summe 


G(x, y) = 2 g,(x, 9) (6. 3) 


ist eine Gitterfunktion, welche wiederum die homogene Gleichung /(#) = 0 be- 
friedigt. In jedem Randpunkt Q, nimmt G(x, y) den Wert 1 an. Denken wir 
uns in (6. 2) g,(x, y) durch G(x, y) ersetzt, so folgt analog fiir jeden innern 
Punkt von G’ 

OG (X,Y) esi. (6. 4) 


Im Spezialfall, wo in (1. 1) der Koeffizient c verschwindet, wo also wegen 
(2. 2) und (2. 3) B = 4 ist, haben wir sogar G(x, y) = 1. 

7. Mit dem Gewonnenen sind wir imstande, zwei Hilfssatze herzuleiten, von 
denen wir in Abschnitt 18 Gebrauch machen werden. 

Hilfssatz 6. Es sei u(x, y) = u(P) eine in G’ definierte Gitterfunktion, 
welche in den Randpunkten Q, (v= 1, 2,..., n), die vorgegebenen Werte U, an- 
nimmt und welche im Innern der homogenen Differenzengleichung l(u) = O geniigt. 
Dann gilt 


u(P) = 3 6 (P) UF (75 


Beweis. Nach Hilfssatz 3 existiert nur eine Gitterfunktion u(P), welche 
gleichzeitig J(u) = 0 befriedigt und am Rande die Werte U, annimmt. Nun 
erfiillt aber die Summe rechter Hand in (7. 1) diese beiden Bedingungen und 
muss daher mit u(P) identisch sein, womit (7. 1) bewiesen ist. 

In Abschnitt 3 wurde jedem Randpunkt Q von G’ ein gewisser Punkt S zu- 
geordnet, der nicht notwendig ein Gitterpunkt ist, der aber doch auf einer 
Gitterseite liegt. Wir erweitern nun den Definitionsbereich von g,(P), indem 
wir analog wie in (3. 2) festsetzen: 


g,(S) = Ay g,(P’) a8 dy £,(P) : ihe /4 


(7. 1) gilt offenbar auch dann noch, wenn man P durch S ersetzt. 
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Hilfssatz 7. Gegeben seien n Zahlen y, mit 0 <n, <1 (yn, = 1). Q bezerchne 


einen beliebigen Randpunkt von G' und S den zugeordneten inneren Punkt. Dann 


existiert eine von der Wahl von Q unabhingige Zahl O mit 


0xsO<1 (73) 
derart, dass 


¥' 7,85) $0 (7.4) 


gilt. 

Beweis. Wenigstens ein 7, ist nach Annahme kleiner als 1. Es sei dies etwa 7). 
Aus (6.4) und (6.1) folgt dann fiir einen beliebigen Gitterpunkt im Innern 
von G’ 


Dy Ny &(P) Sm &(P) + &a(P) + +++ + Bn(P) = G(P) + 8x(P) (m1 — 2) 
f = 1-—(1—m) g(P) <1. 


Im Falle, wo P ein Randpunkt von G’ ist, kénnen wir dieselbe Abschatzung 
durchfiihren. Nur muss in der letzten Ungleichung das Gleichheitszeichen zu- 
gelassen werden, da ja g,(P) = 0 sein kénnte. Wenden wir die Ungleichung auf 
P' und P" an und beachten, dass nach Abschnitt 3 wenigstens einer von ihnen 
ein innerer Gitterpunkt ist, so folgt aus (7. 2) wegen A, => 0,4, 20,4, +4,=1 


sofort » 7, g,(S) <1. Da nur endlich viele Randpunkte Q mit je einem zuge- 
v=1 


ordneten Punkt S vorhanden sind, so existiert die im Hilfssatz 7 angegebene 
Schranke 9 < 1. 

8. In diesem Abschnitt rechtfertigen wir die Voraussetzung, wonach die 
Randkurve C von Ecken (und Spitzen) frei sein soll. 

Die Schwierigkeiten, welche Ecken bieten, rithren nicht von der Randbe- 
dingung (1. 2) her. Denn um eine samt ihren partiellen Ableitungen stetige 
Lésung v(x, y) zu erméglichen, geniigt es, den Koeffizienten k(s) und A(s) in 
allfalligen Eckpunkten passende Unstetigkeiten beizulegen. 

Auch wird der Rechner vor solchen Singularitaéten der Randkurve kaum 
zurtickschrecken. Es wird ihm ohne Miihe gelingen, geeignete Gitterrandpunkte 
und brauchbaren Ersatz fiir dahinfallende Randbedingungen in der Nahe von 
Ecken zu finden. 

Vielmehr sehen wir uns deshalb gezwungen, Ecken in der Randkurve auszu- 
schliessen, weil unsere Untersuchung der Approximation einer Gitterlésung an 
die exakte Losung vor noch nicht tiberwundene Schwierigkeiten gestellt wiirde. 
An zwei Beispielen der Figur 2 weisen wir auf die Besonderheiten hin, die bei 
Ecken vorliegen k6nnen: 

In Figur 2a kénnen wir vom Gitterrandpunkt Q kein Lot auf die Randkurve 
fallen und daher eine Randbedingung zweiter oder dritter Art nicht approxi- 


~ 
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mieren. In Figur 2b ist eine Ecke mit kleinem Winkel « gezeichnet. Hier 
chaufen » sich die Gitterrandpunkte derart stark, dass nicht einzusehen ist, wie 
die erforderlichen Differenzengleichungen aufgestellt werden kénnten. 


ak 
+— + 
C 
“| $Q_ | 
a 
oeaeice 
a r) 
Fig. 2 
Die Schwierigkeiten der Approximation, wenn die Randkurve C Ecken besitzt. 


Wahrend wir uns bei stetiger Tangente die eindeutige Zuordnung der Git- 
terrandpunkte zu den Punkten von C dadurch sichern konnten, dass wir die 
Gitterseite 4 hinreichend klein voraussetzten, versagt, dieser Weg bei Ecken 
vollkommen. Denn wie fein wir auch ein Gitter wihlen, so werden ber Ecken 
keinerler Erleichterungen eintreten. 

Nur im Falle eines Stufenrandes, wo die Randkurve C stiickweise parallel 
- zu den Gitterseiten verlauft, bilden Eckpunkte fiir unsere Betrachtungen kein 
Hindernis. 


§ 2. Der beim Ubergang zu den Differenzengleichungen 
begangene Fehler 


9. Bei der im folgenden durchzufiihrenden Fehlerabschatzung wird es sich 

darum handeln, den Unterschied der Lésung v(x, y) der Differentialgleichung 
(1. 1) mit der Randbedingung (1. 2) gegentiber der Losung u(x, y) des zugeord- 
neten Systems der Differenzengleichungen (4.1) zu bestimmen. Fiir diese 


Differenz setzen wir 
(x, y) = (4%, y) — u(x, y) - (9-1) 


6(x, y) kann natiirlich nur in den Gitterpunkten von G’ gebildet werden. Es 
ist zu erwarten, dass sich 6(x, y) beim Ubergang zu einem feineren Gitter ver- 
kleinert. Es wird unser wichtigstes Anliegen sein, zu priifen, wie | d(x, y)| von 
der nach Null strebenden Gitterseite h abhangt. 

Vorerst bestimmen wir die Restglieder, welche beim Ubergang von der Dif- 
ferentialgleichung zu den Differenzengleichungen auftreten. Zu diesem Zweck 


ZAMP III/12 
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sei w(x, y) eine beliebige Funktion, die in G samt ihren partiellen Ableitungen 
bis zur vierten Ordnung stetig ist. Wie in Figur 1 greifen wir einen innern 
Punkt P, und seine vier Nachbarpunkte im Gitter heraus und bezeichnen die 
Werte von w(x, y) und der Ableitungen in P, kurz mit w, 0w/dx , ..., die Werte 
von w(x, y) in den Nachbarpunkten aber mit w,, ..., @y. Dann haben wir 


W, | ; Ow iP MOET, he 08 ht (% 

ie ee oe T Ol Oa® 31 Oa + 41 LY 

Ws | 4a: 0 3 (9 2) 
| \ ow ht. S0%e AY One eee | 

ws ee Oye OT Gye ay age srs Lv 


i ance nagmmmagtae, 


wo Y,, Y% und ¥,, Y, gewisse mit P, variable Zwischenwerte der vierten par- — 


tiellen Ableitung von w(x, y) nach x bzw. nach y bezeichnen. Bricht man in 
(9. 2) schon bei Gliedern dritter Ordnung ab, so wollen wir die entsprechenden 
Zwischenwerte der partiellen Ableitungen dritter Ordnung mit ®,, ®, bzw. 
mit ®,, ®, bezeichnen. 

Durch passende Kombination der Relationen (9. 2) folgt 


0?2w 02w 1 h2 


‘apts ) oye An (W, + Wy + W, + Wy —4v) — =p CPT shee 
und 

Ow W,— Ws I OW = Wa — W, h? 

Oe ah ag Pa Oly | a gee 


Fiihren wir noch zur Abkiirzung fiir den Differentialausdruck in (1. 1) die 
Bezeichnung 
Top 0?w 0?w Ow = 


l(w) = pet cage aye + cw (9.3) 
ein, so folgt 
1 nee 
lw) = he Uw) — 3 zz (Py +--+) + a (®, + ®,) + 6 (D+ ®,)| » (9. 4) 


wo /(w) den in (4. 1) definierten linearen Ausdruck bezeichnet. Das Restglied | 
ohne den Faktor h® bezeichnen wir mit 


¥ rad 
Ry = ap [a it $F) +4 (+O) +5, 4+H)]. 5) 


R,, ist ein lineares Funktional von w(x, y). Nunmehr kann (9. 4) in der Gestalt 


h? l(w) = l(w) + h* R, (9. 6) 
geschrieben werden!), 


1) Mit anderen Bezeichnungen hat S. Gerscucorrn [10], S. 375, die Relation (9. 6) zusammen 
mit dem Restglied angegeben. 


i 


7 


{ 
> 
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Wenden wir (9. 6) auf die Lésung v(x, y) von (1. 1) an, was wegen der in 
Abschnitt 1 gemachten Annahmen zulassig ist, so folgt 


h?g=l(v) + h*R, 


oder nach Subtraktion der fiir u(x, y) giiltigen Relation (4. 1) bzw. (2. 2) wegen 
(9. 1) 
(6) = 1(v—u) = —h*R,. (Oren) 


Diese Relation gilt fiir jeden innern Gitterpunkt von G’. Der Faktor R, ist von 
der Gitterseite / indirekt abhangig. Da aber alle in R,, eingehenden Koeffizien- 
ten und Ableitungen beschrankt sind, so existiert eine von # unabhangige, posi- 
tive Schranke R*, so dass 

|1(5)| <h4 Re (9. 8) 
gesetzt werden darf. 

10. Etwas mehr Miihe bereitet die Diskussion der Fehler, die durch Appro- 
ximation des Randes C durch den Streckenzug C’ zustande kommen. In Figur 3, 
welche einen kleinen Ausschnitt von Figur 1 wiedergibt, betrachten wir einen 
beliebigen Gitterrandpunkt Q auf C’. Die durch Q laufende Normale von C 
ordnet Q einen Punkt Q auf C zu sowie auch einen Punkt S auf der nachsten 
von der Normalen geschnittenen 
Gitterseite im Innern von G’. Wie 
in Abschnitt 3 seien die Endpunkte 
dieser Gitterseite mit P’ und P” be- 
zeichnet. Setzen wir QQ = 0 h, so gilt 
wegen der in Abschnitt 2 dargelegten 
Konstruktion von C’ 


(10. 1) 


wobei 0 positiv zu nehmen ist, wenn 
Q innerhalb von C liegt, und negativ 
im umgekehrten Fall. Uber QS = oh 
liegt bereits die Ungleichung (3. 1) 
vor. Ferner wird mit (3. 2) das Teil- 
verhaltnis des Punktes S auf P’P" 
festgelegt, so dass P’S =/,h und 
P'S: =}, h za setzen ist. 

Es sei w(x, y) wie in Abschnitt 9 
eine stetige Funktion in G, deren Ableitungen passenden Beschrankungen 
unterworfen sind. Die Werte, welche w(x, y) in den Punkten Q, Q und S 
annimmt, seien der Reihe nach mit W, W und @ bezeichnet. Um mit w(x, y) 
langs dem Gitterrand C’ die in (4. 1) vorkommende Linearform 


Fig. 3 


Approximation der Normalableitung. 
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L(W)=Kw—-HW (10. 2) 


bilden zu kénnen, haben wir zunachst zu beachten, dass w nicht mit @ tiber- — 


ee arte 


einstimmt, trotzdem sich beide Werte auf den Punkt S beziehen. Denn w muss ~ 


nach Vereinbarung von Abschnitt 3 durch lineare Interpolation aus w(P’) und 
w(P") bestimmt werden. Die Differenz @ — w ist folglich der Interpolations- 
fehler, der im Falle, dass P’P” parallel zur x-Achse verlauft, lautet 


2 (%¥,— Xp) (¥;— *pr) / 2h’ 
ie ae ( Pe 


2! Ox? 


Beriicksichtigen wir auch eine Interpolation in der y-Richtung, so folgt allgemein ~ 


D—-w=h* J, (10. 3) 


wo /,, eine variable Grésse ist, fiir die 


| 02w | 


i! 
bie ys liscrp 


Oy? | 


gilt. Diese obere, von h unabhangige Schranke sei mit /* bezeichnet. 
11. Um die in (1. 2) vorkommende Linearform 


die jetzt fiir den Randwert W im Punkte Q geschrieben ist, mit L(W) in Bezie- 
hung bringen zu kénnen, driicken wir W und w mit einer Taylor-Entwicklung 
durch die Werte und Ableitungen am Rande aus: 


wo 
T= 0O0+0 
gesetzt ist und wo 7%, 7, gewisse Mittelwerte der Normalableitung zweiter Ord- 


nung und y; einen Mittelwert erster Ordnung bezeichnen. 


Fiir Punkte auf C, gewinnen wir wegen k(s) = 1 zunadchst aus (10. 2), 
(10. 3), (11. 2) und (11, 3) 


LW) = (1-H) W+ @h—- Hoh) So 4m |S y— He x — J]. 


on 


Hierin ersetzen wir H durch den in (3.8) stehenden Ausdruck und 0W/dn 


gemiss (11.1) durch L(W) + h(s) W. Ferner fiihren wir zur Abkiirzung an 


a ) (10. 4) | 


—h(s) W , (11.1) 9 
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Stelle der eckigen Klammer die Grésse T,, ein. Dann folet wegen t =0 +0 


 L(W) = oh L(W) — h? [9c h(s) L(W) + 9c h%\(s) W — Ty] auf Ci. (11. 4) 
Auf C3 aber schliessen wir aus (10. 2) und (11. 2) wegen K = 0 und H=oh 
L(W) =—chW=-—chW — och? 4, 
und, wenn wir noch (11. 1) mit k(s) = 0, h(s) = 1 beiziehen, 
L(W) =ohL(W)— och? yx, aufCs. (11. 5) 


Damit ist in beiden Fallen L(W) durch L(W) ausgedriickt. 

12. Wenden wir (11. 4) und (11.5) auf die spezielle Funktion v(x, y) an, 
welche die Lésung der Gleichung (1. 1) mit den Randbedingungen (1. 2) be- 
zeichnet, so erhalten wir 


L(V) =ohf\(s) —h? [oo h(s) f(s) + eo hs) V— T,] auf Ci, 
L(V) =chf(s) fh? oo x, aut C2. 


Subtrahieren wir davon die in (3.5) fiir u(x, y) hergeleitete Randbedingung 
L(U) =F=ohf{i(s), so erhalten wir wegen (9. 1) 


L(6) = —h® {9 o h(s) f(s) + 9.0 h%(s) V — T,] auf Ci, (12, 1) 
L(6) = —h? oo x5 ah Ce (1252) 


Fiir die absoluten Betrage der mit * multiplizierten Ausdriicke existiert, ahn- 
lich wie in (9. 8), eine von # unabhangige, gemeinsame obere Schranke S¥, so 
dass wir 

LE (6) |b? Sk waut.Co (12333) 
schreiben diirfen. 

13. Nach diesen Vorbereitungen gelingt es uns nun, das von S. GERSCH- 
GORIN [10] fiir den Spezialfall der Randbedingungen erster Art geschaffene 
Majorantenverfahren auf Probleme mit anderen Randbedingungen zu verall- 
gemeinern. 

Uber die Fehlerfunktion 6(x, y), die in (9. 1) definiert wurde, haben wir bis 
jetzt die Ungleichungen (9. 8) und (12. 3) hergeleitet. Um von hier aus zu einer 
Abschatzung von | d(x, y)| zu gelangen, fiihren wir die Hilfsfunktion 2(x, y) 
als Losung des folgenden Randwertproblems ein: 


> 022 O22 Oz Oz . 

i(z).= 4 +a l= -¢f=—1> mG, | 
LE ag y (13. 1) 

L(z) = As) oe —h(s)z=—1 aufC. | 
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Nach den in Abschnitt 1 getroffenen Voraussetzungen existiert eine eindeutig 
bestimmte Lisung z(x, y), die samt ihren partiellen Ableitungen bis zur vierten 
Ordnung in G stetig ist). ; 
Wir versuchen nun, | 6(x, y)| durch eine von 2(x, y) abgeleitete Hilfsfunk- 
tion 
w(x, vy) = m2(x, ¥) (13. 2) 
nach oben abzuschatzen. Hierin bezeichnet m einen passend zu bestimmenden, 


von x, y unabhangigen, positiven Parameter. Wir stiitzen uns dabei auf den 
Hilfssatz 5, wonach 


| (x, y) | S wl, 9) (13. 3) 

bereits aus den Ungleichungen 
Iw) <—|1(6)| im Innern von G’, (13. 4) 
L(W) =—|L(6)| auf C’ (13.3 


folgt. Es wird also unsere Aufgabe sein, durch passende Wahl von m die Un- 
gleichungen (13. 4) und (13. 5) zu erzwingen. 

Um (13. 4) zu erfiillen, greifen wir auf die Formel (9. 6) zuriick, welche auf 
die in (13. 2) definierte Funktion w(x, y) wegen der 2(x, y) und ihren partiellen 
Ableitungen zukommenden Stetigkeitseigenschaften anwendbar ist. Wir er- 
halten 

I(w) = h? l(w) — h* R, = mh? l(z) — mh* R, 
oder wegen (13. 1) 
lw) = —mh? (1+ h? R,). 


Nun sei R¥ eine obere Schranke fiir | R,| im Gebiet G und h bereits so 


klein angenommen, dass 


h? R* < 


nr] 


ist. Dann folgt 


Uw) < — 


nl 


mh*, (ise 
Um also (13. 4) zu erhalten, geniigt es, wegen (9. 8) 


i mh? < — h* R* 
oder 
me Zh* Ke (13.77) 
zu setzen. 


1) In der Definition von z weichen wir von S. GERSCHGORIN ab, indem wir die Randbedingung 


inhomogen ansetzen. Dadurch ersparen wir uns eine miihsame und uniibersichtliche Abschatzung des 
in (13. 2) eingeftihrten Faktors m. 
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14. Zur Herleitung von (13. 5) bentitzen wir (11. +k) (1 Seb \tundter- 
halten 


LW) =—mho{i—hlens)— eh) Z+ 47] auf C!, 
LW) =—mho{1+ he x,(z) aut Cy. 


h sel bereits so klein gewahlt, dass die absoluten Betrage der beiden ge- 
schweiften Klammern tiberall auf Ci, bzw. Ci mindestens den Wert 1 /2 be- 
sitzen. Dann folgt wegen o = 1 in beiden Fallen gemeinsam 


LW) < —+ mh. 
Damit (13. 5) erfiillt ist, haben wir wegen (12. 3) bloss 


1 
7 mh < —h? S* 


oder 
mh 2 2h os (14. 1) 
zu setzen. 
Mit (13. 7) zusammen miissen wir somit 
teen Max (2h R255) (14. 2) 


fordern. Fir hinreichend kleines / existiert somit eine von / unabhangige, posi- 
tive Konstante M, so dass zur Erfiillung von (13. 4), (13.5) und damit von 
HL3s 3) 

m=h M (14.3) 
gesetzt werden darf. 

Fihren wir noch zur Abkiirzung die von der Feinheit des Gitters unabhan- 
gige Funktion ¢(x, y) = M z(x, y) ein, so kénnen wir das in (13. 3) erhaltene 
Ergebnis in dem folgenden Satz formulieren: 

Satz 1. List man an Stelle der gegebenen Differentialgleichung (1. 1) mit den 
_ Randbedingungen (1.2) das System der Differenzengleichungen (4.1) fiir ein 
quadratisches Gitter mit der Seite h auf, so gilt fiir den Unterschied 0(x, y) der 
beiden Losungen in den Gitterpunkten die Ungleichung 


| d(x, ¥) | Sht(x, y), (14. 4) 


wo t(x, y) eine nur vom gestellten Problem abhdangige, in G stetige Funktion 
bedeutet. 

Geht man demnach zu immer feinern Gittern tiber, so strebt der Fehler 6(x, v) 
mindestens linear mit h gegen Null, und die Gitterlosung konvergiert gleichmdssig 
gegen die exakte Lisung des gestellten Problems. 

15. Bei einem praktischen Beispiel wird man selten in der Lage sein, ¢(x, y) 
zuverlassig genug zu bestimmen. Uberdies diirfte sich der Aufwand kaum lohnen. 
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Es ist jedoch nicht unwichtig, zu wissen, von welchen Parametern und Funk- 
tionen ¢(x, y) abhangt, um in Fallen schlechter Konvergenz die besten Mass- 
nahmen zur Beschleunigung der numerischen Rechnungen treffen zu kénnen. 

Sollte zum Beispiel eine guadratische Konvergenz erwiinscht sein, das heisst 
in (14. 4) der Faktor A? an Stelle von h stehen, so ersieht man aus (14. 2), dass 
S* linear mit h gegen Null streben miisste. Nun ist S% eine gemeinsame 
obere Schranke fiir die Betrage der mit h? multiplizierten Ausdriicke in (12. 1) 
und (12. 2). Die erste dieser Relationen steht fiir eine Randbedingung dritter, 
die zweite fiir eine solche erster Art. In beiden Fallen miisste man also die 
Approximation an die exakte Randbedingung (1. 2) verbessern. Dies ist auf 
Kosten einer vermehrten Rechenarbeit ohne weiteres méglich. Fiir eine Rand- 
bedingung erster Art empfiehlt sich die von L. CoLLatz[1] vorgeschlagene 
Interpolationsmethode. Um das Ziel auch im Fall einer Randbedingung zweiter 
oder dritter Art zu erreichen, miisste man den in (3. 3) eingefiihrten Differen- 
zenquotienten durch einen finiten Ausdruck zweiter Ordnung ersetzen. Geo- 
metrisch bedeutet dies: Man miisste die Ableitung am Rande durch eine Drei- 
punkteformel der numerischen Differentiation approximieren?). 

Unser Ergebnis zeigt auch, dass die gelegentlich gedusserte Meinung, man 
miisse vor allem den Approximationsgrad der Differenzengleichungen (2. 3) im 
Innern des Gebietes erhdhen, jedenfalls dann nicht zu Recht besteht, wenn es 
darum geht, den Abbrechfehler (truncation error) zu vermindern. 


§ 3. Das numerische Verfahren zur angeniherten Auflésung 
des zugeordneten Systems von Differenzengleichungen 


16. Das System von Differenzengleichungen (4.1) enthalt ebenso viele 
Gleichungen, wie das gewahlte Gittergebiet G’, einschliesslich des Randes C’, 
Punkte besitzt. Diese Zahl ist in praktischen Beispielen ziemlich gross, und es 
ist daher keine Rede davon, das System algebraisch auflésen zu kénnen. Es 
wird nichts anderes iibrigbleiben, als nach einem geeigneten Approximations- 
verfahren zu suchen. 

In dem Spezialfall, wo iiberall auf C’ eine Randbedingung erster Art vor- 
hegt, fallen die auf den Rand beziiglichen Differenzengleichungen in (4. 1) fort, 
da ja die Randwerte fest sind. Fiir diese spezielle Aufgabe empfiehlt sich das 
vor einigen Jahren von R. V, SourHweELt [21] ersonnene Relaxationsverfahren, 
dessen Uberlegenheit iiber frithere Methoden mehrfach bestitigt wurde, vgl. 
[6], [7], [81, (221, [23}2). , 


1) Die von L, Fox [8] vorgeschlagene Dreipunkteformel eignet sich anscheinend nur fiir einen 
rechteckigen oder einen Stufenrand. 

*) Wir sehen hier und im folgenden von der Verwendung von Rechenautomaten ab. Selbstver- 
standlich ergeben sich fir das Maschinenrechnen etwas andere Gesichtspunkte. Wir verweisen 
hiezu auf die kiirzlich von S. P. Franxex [9] und von E, Stiret [22] erschienenen Arbeiten. 
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Zur numerischen Aufldsung von (4. 1) fiir den allgemeineren Fall, wo auch 
am Rande Differenzengleichungen auftreten, liegt es nun nahe, das Relaxa- 
tionsverfahren von SOUTHWELL in geeigneter Weise zu erginzen. Dies ist der 
Grundgedanke unseres Vorgehens in Abschnitt 17. 

Das Relaxationsverfahren geht wie jedes Iterationsverfahren von Roh- 
werten w = w(x, y) aus, welche hier ganz beliebig sein diirfen. Die Differenzen- 
gleichungen (4. 1) sind fiir die Rohwerte natiirlich nicht erfiillt. Wir setzen aber 


iw) —y =71, (16. 1) 
LW)-F=R. (16. 2) 


y = r(x, y) ist eine Gitterfunktion, die fiir alle Punkte im Innern von G’ defi- 
niert ist. R = R(x, y) ist nur fiir Randpunkte erklart. Die Zahlen ry und R 
heissen die Restduen von w in den Gitterpunkten. Die Residuen der Funktion 
u = u(x, y), welche das Gleichungssystem (4. 1) lost, sind offenbar Null. Umge- 
kehrt kann man sich aus dem Wert der Residuen einer Versuchsfunktion 
w = w(x, y) ein ungefahres Bild von dem «Fehler» w — u machen. 

Beim Randwertproblem erster Avt kOnnen die Residuen R am Rande von 
vorneherein als Null vorausgesetzt werden. Es handelt sich dann darum, die 
Absolutwerte der Residuen 7 in allen inneren Punkten schrittweise kleiner und 
kleiner zu machen. Hiefiir wurde eine Reihe konvergenzbeschleunigender 
Kunstgriffe geschaffen, die in der Literatur ausfiihrlich beschrieben sind (vgl. 
[9], [21], [22], [23]). Die Reduktion der Residuen wird kurzweg Relaxation ge- 
nannt. Dass dabei w(x, y) gegen u(x, y) konvergiert, ist unter verschiedenen 
Voraussetzungen, bald iiber die Differentialgleichung, bald iiber die Reihen- 
folge des Vorgehens, mehrfach bewiesen worden (siehe [3], [4], [11], [12], [15], 
[18], [22], [23], [24]). 

17. Treten nun Randbedingungen zweiter oder dritter Art hinzu, so hat man 
ausser den Residuen 7 auch die Residuen R in Randpunkten zu reduzieren. ° 
Nach (4. 1) hangen 7 und R sowohl von den Funktionswerten w(x, y) in innern 
Punkten als auch von denen in Randpunkten zugleich ab. Man kann somit die 
Residuen am Rande nicht unabhangig von den Residuen im Innern gegen Null 
streben lassen. 

Es liegt nahe, in der Weise vorzugehen, dass man beim Abandern von 
~ w(x, y) in randnahen Punkten gleichzeitig die Randwerte so mitverandert, dass 
R.= 0 erreicht wird. Der Rechner diirfte so am raschesten zum Ziele gelangen. 
Um aber die Konvergenz des Verfahrens leichter beweisen zu kénnen, empfiehlt 
sich fiir uns eher ein getrenntes Vorgehen, das sich auf den folgenden beiden 
Schritten aufbaut: 

1. Schritt. Ausgehend von irgendeiner Gitterfunktion w(x, y) in G’ bringe 
man die Residuen (16. 1) in allen inneren Punkten nach dem iiblichen Relaxa- 
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tionsverfahren zum Verschwinden, lasse aber die Randwerte W(x, y) unver- 
andert?). 

2. Schritt. Durch Anpassung der Randwerte W(x, y) mache man die Resi- 
duen (16. 2) in allen Randpunkten zu Null, 

Beim ersten Schritt werden die Residuen im Innern zwar Null. Am Rande 
treten aber von Null verschiedene Residuen auf. Das Umgekehrte erfolgt beim 
zweiten Schritt. Daher hat man die beiden Schritte abwechslungsweise fort- 
zusetzen. Im folgenden Abschnitt wird gezeigt, dass dieses alternierende Ver- 
fahren konvergent ist. 

18. Satz 2. Bringt man, ausgehend von einer beliebigen Gitterfunktion w(x, y), 
abwechslungsweise die Residuen (16.1) mit Hilfe des Relaxationsverfahrens und 
die Residuen (16. 2) durch geeignete Abanderung der Randwerte zum V erschwin- 
den, so konvergiert die Gitterfunktion gegen die gesuchte Lésung u(x, y) des Glet- 
chungssystems (4. 1). 

Beweis. Beim ersten Schritt gehen beliebig angenommene Rohwerte w(x, v) 
mit den Randwerten W in eine neue Gitterfunktion w(x, y) tiber, deren Rand- 
werte W, mit W iibereinstimmen. Nur ist jetzt y = 0 in (16. 1), das heisst es gilt 


Iw) =y- 
Fiir die Residuen am Rande haben wir zu setzen 
R=L(W,)-F=Kw,—HwW-F. (18. 1) 


Beim zweiten Schritt bringen wir die Residuen am Rande zum Verschwinden, 
indem wir zu den Randwerten W passende Gréssen AW addieren. Die so ge- 
wonnene neue Gitterfunktion sei mit w,(*, v) bezeichnet. Sie stimmt in allen 
inneren Punkten mit w,(*, vy) iiberein, so dass wir 


0=L(W,) -F=Kuv,—-HW,-F=Kw,—H(W+AW)-—F (18. 2) 
haben. Ein Vergleich mit (18. 1) fiihrt zu 
Rai AW, (18. 3) 


Fir Punkte auf Cj sind hierin AW und R gleich Null zu setzen. 
Beim letzten Vorgang sind im Innern neue Residuen entstanden, welche bei 
nochmaliger Anwendung des ersten Schrittes zwar verschwinden, jedoch ge- 


wisse Anderungen Aw = Aw(x, y) an w (x, y) in den inneren Gitterpunkten zur 
Folge haben. Wir setzen daher 


W(x, y) = w(x, y) + dw=w,(x, y) + dv, Wa=We=W+ AW. 


1) Praktisch muss man wegen der beschrainkten Anzahl von Dezimalen die Reduktion der 
Residuen abbrechen, sobald sie, in der Einheit der letzten mitgefiihrten Dezimalen ausgedriickt, 
kleiner als 2 sind und keine Blécke von vorwiegend positiven oder negativen Residuen vorhanden 
sind. Man vergleiche hiezu die Betrachtungen von E, Stirret [22]. 


a a a Te ae ee ty A lr TT 
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Die Anderungen Aw lassen sich mit Hilfe der in Abschnitt 6 eingefiihrten 
Greenschen Zahlen folgendermassen berechnen: 

Da w, mit den Randwerten W die Gleichungen /(w) = y erfiillt und dasselbe 
auch von wz mit den Randwerten W + AW gilt, so muss die Differenz der bei- 
den Gitterfunktionen eine neue Funktion sein, welche der homogenen Gleichung 
L(ws — wW,) = 0 geniigt. Ihre Werte im Innern sind gleich Aw, ihre Werte am 
Rande gleich AW. Bezeichnet man die Gitterpunkte von G’ mit Q,, ..., Q,, und 
ordnet ihnen die Randwerte AW,, ..., AW,, zu, so folgt aus Hilfssatz 6 fiir jeden 
innern Gitterpunkt P 


Aw(P -» AW, g(P (18. 4) 


Schreibt man noch (18. 3) mit dem Index », so geht (18. 4) wegen H, > 0 in 


Aw(P) = S7 a g,(P) (18. 5) 


uber, wo fiir Punkte, die auf Cj liegen, R, = 0 zu setzen ist. 

Weiter berechnen wir nach Formel (18. 1) die neuen Residuen k’, die beim 
Ubergang von w, zu w, + Aw = w,am Rande entstehen miissen. Beachten wir, 
dass w, fiir alle inneren Punkte gleich w, und am Rande gleich W + AW ist, 
so folgt 

Rave) =F = we Aw) AW PAW) =F 


oder wegen (18.1) und (18. 3) 
Ko= K Aw. (18. 6) 


Fiir Randpunkte auf Cj verschwinden K und R’. Fiir Randpunkte auf Cj ist 
K = 1und Aw ein Betrag, der nach den Ausfiihrungen von Abschnitt 3 folgen- 
dermassen zu berechnen ist: Jedem Randpunkt Q,, von Cj ist ein Punkt S,, 
zugeordnet, der nicht notwendig ein Gitterpunkt zu sein braucht. Entspre- 
chend ist der Faktor Aw = Aw(S,) in (18. 6) nach (3. 2) durch lineare Inter- 
polation aus den Werten von Aw in P’ und P" zu bestimmen. Beniitzen wir 
(18.5) mit P= S,, (was nach der dem Hilfssatz 6 folgenden Bemerkung zu- 
lassig ist), so besagt (18. 6) fiir einen Punkt Q,, von Cy 


R’, = Av(S = SH a ~ g,( (18. 7) 


Mit R, und Rj bezeichnen wir dasjenige unter den Residuen R, bzw. Kk, 
mit maximalem absolutem Betrag. Dann folgt aus (18. 7), wegen (3. 8) und 


(6. 1), 
ig ea g, | S) S| Rel 3 ae 61S) (18. 8) 
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wo 6, gleich 1 oder 0 zu setzen ist, je nachdem ob Q, zu C; oder zu Cy gehort. 
Setzen wir noch 

ae =), (y= 1,2, <2., 8) 
so gilt wegen (3. 8) und (3. 9) 0 <7, <1 (yn, + 1). Daher lasst sich Hilfssatz 7 
auf die Summe rechter Hand in (18. 8) anwenden, so dass wir 


|R,| SO|R, 


(O <1) (18.9) 


erhalten. Dies bedeutet: Das Maximum des absoluten Betrages der Residuen 
am Rande ist kleiner geworden. Nach N-maliger Anwendung der beiden 
Schritte hat es sich mindestens um den Faktor OY verkleinert. Daher streben 
die Residuen am Rande gegen Null. 

Dass gleichzeitig auch die Residuen im Innern gegen Null streben, schliessen 
wir aus (18. 3). Danach miissen die Anderungen AW proportional mit den Resi- 
duen am Rande gegen Null konvergieren. Diese AW sind es aber, welche gemass 
(16. 1) in randnahen inneren Punkten gewisse Residuen erzeugen, die wegen 
der Beschrankungen, denen die Koeffizienten «; in (16. 1) unterliegen, ebenfalls 
eine Nullfolge bilden miissen. 

Schhiesslich folgt aus der Eindeutigkeit und der Stetigkeit der Lésung von 
(4. 1), dass die Naherungsfunktion w(x, y) gegen die exakte Lésung u(x, y) 
konvergiert, womit der Satz 2 bewiesen ist. 

19. Die im Beweis von Satz 2 gefundenen Relationen gestatten eine Dis- 
kussion dariiber, mit welcher Geschwindigkeit das Approximationsverfahren 
zur numerischen Auflésung des Gleichungssystems (4. 1) konvergiert. Offenbar 
ist die Geschwindigkeit um so grésser, je kleiner der in (18. 9) stehende Faktor 0 
ist. Die Betrachtung von (18. 8) fiihrt zu den folgenden Fallen, in denen @ ver- 
haltnismassig klein ausfallt und in denen folglich die Geschwindigkeit der Kon- 
vergenz befriedigend ist: 

a) Cy enthalt einen nicht zu kleinen Anteil an Randpunkten von G’, das 
heisst, es existieren gentigend viele Randpunkte mit einer Randbedingung 
erster Art. Ftir solche Punkte verschwindet 6, in (18. 8) und 9, in (7. 4), so dass 
die Behauptung tiber 0 unmittelbar folgt. 

b) Die in (6. 3) definierte Summe G(x, y) ist fiir innere Gitterpunkte wesent- 
lich kleiner als 1. In diesem Falle folgt die Kleinheit von @ ohne weiteres aus 
dem Beweis von Hilfssatz 7. Damit aber G(x, y) hinreichend klein ausfallt, 
gentigt es nach den Betrachtungen in Abschnitt 6, dass die in (2. 3) eingefiihrte 
Grésse B = 4 — ch? merklich grésser als 4 sei. Somit muss der Koeffizient c in 
der Differentialgleichung (1.1) wesentlich kleiner als Null vorausgesetzt werden. 
Man erkennt jedoch, dass der Vorteil bei kleiner Gitterseite # verlorengeht. 

c) Der in (3. 6) eingefithrte Koeffizient H ist in den Punkten von Cj er- 
heblich grésser als 1. Wegen (3. 8) ist H = 1 + o h h(s). In diesem Fall muss also 
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der Koeffizient 4(s) der Randbedingung (1. 2) auf C, wesentlich grosser als Null 
sein. Wie bei b fallt aber der Vorteil der raschen Konvergenz dahin, wenn die 
Gitterseite # hinreichend klein ist. 
Zusammenfassend kann man sagen: Bei grobem Gitter kann die Geschwin- 
digkeit der Konvergenz aus mehreren Griinden gut sein. Ist man aber zur Erhohung 
der Genauigkeit auf ein feineres Gitter angewiesen, so sinkt die Geschwindigheit 
der Konvergenz mit abnehmender Gitterseite h.‘Sie bleibt befriedigend, wenn der 
Anteil der Randpunkte mit einer Randbedingung erster Art hinreichend gross ist. 

Bei einem konkret vorliegenden Beispiel lasst sich allerdings der Konver- 
genzgrad durch geschickte Massnahmen erheblich steigern. Man verbessert im 
_ Laufe des Verfahrens die Randwerte so, dass die Residuen R, am Rande nach 
der Reduktion der Residuen im Innern nicht einseitig positiv oder negativ aus- 
fallen. Dann lasst sich die Abschatzung in (18. 8) verscharfen, so dass sich ein 
_ kleinerer Wert fiir @ ergibt. 

20. Bis jetzt haben wir in § 3 mit einem festen Gitter gearbeitet. Bei einem 
numerischen Beispiel wird man aber die Gitterseite schrittweise verkleinern. Ein 
_ grobes Gitter mit bloss etwa 30 bis 60 Gitterpunkten dient dazu, dem Rechner 
eine erste Orientierung tiber den Verlauf der Losungsfunktion zu geben. Bei 
der Approximation der Randkurve sind die Regeln von Abschnitt 2 und 3, 
 soweit es die beschrankte Punktzahl zulasst, einzuhalten. Zur Erleichterung 
wird man anfanglich die Gitterfunktion bloss zwei- oder dreistellig wahlen. 
Die Residuen sind zu reduzieren, bis sie nur noch eine oder zwei Einheiten der 
- letzten mitgefithrten Dezimalen betragen und bis nirgends mehr «gefahrliche » 
Blécke mit vorwiegend positiven oder negativen Residuen vorhanden sind?). 

Die Verfeinerung des Gitters wird am besten dadurch vorgenommen, dass. 
man die Gitterseite h halbiert. Auf diese Weise kommen zu den alten Gitter- 
punkten noch ungefahr dreimal so viele neue hinzu. Die Approximation der 
Randkurve hat von neuem zu erfolgen, wobei die Regeln von Abschnitt 2 und 
3 im allgemeinen leichter eingehalten werden kénnen als beim gréberen Gitter. 
- Als Rohwerte fiir die neue Gitterfunktion bentitze man in den alten Gitter- 
- punkten die bei der ersten Approximation erhaltenen Werte. In den neuen 

inneren Gitterpunkten aber ermittle man die Rohwerte aus den schon bekann- 
ten durch Interpolation oder besser durch geschickte Beniitzung der Relation 
I(u) = y. Dann beginnt wieder die Reduktion der Residuen usw. 

Bei fortschreitender Verfeinerung des Gitters muss die Rechengenauigkeit 
durch Hinzunahme weiterer Dezimalstellen erhéht werden, da sonst die Ab- 
rundungsfehler eine verhangnisvolle Rolle spielen werden. In welchem Ausmass. 
dies zu geschehen hat, wird in Abschnitt 22 noch untersucht werden. 
21. Besonders giinstig steht es um die Méglichkeit von Rechenkontrollen. 
-Wahrend der Reduktion der Residuen kénnen sich Rechenfehler bei der be- 
trachtlichen Zahl der in den Rechenprozess eingehenden Werte verhaltnis- 


1) Vgl. Fussnote 1 auf Seite 186. 


| 
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missig leicht einstellen. Zum Zweck ihrer Eliminierung berechne man von Zeit 
zu Zeit die Residuen aus (16.1) und (16. 2) neu. Fast unfehlbar wird dabei 
jedes Versehen entdeckt. f 

Aus diesem Grunde zeichnet sich das Differenzenverfahren und das darauf 
fussende Relaxationsverfahren vor anderen numerischen Verfahren dadurch 
aus, dass man im Laufe des Approximationsprozesses immer wieder von vorne 
beginnt und nicht auf schon efhaltenen, méglicherweise mit Fehlern behaf- 
teten Werten endgiiltig aufbauen muss. 

22. Wie schon in Abschnitt 20 angedeutet wurde, spielen die Abrundungs- 
fehler eine erhebliche Rolle. Da man nur eine beschrankte Zahl von Dezimal- 
stellen mitfiihren kann, so wird jede Gitterlésung u(x, y) notwendig mit ge- 
wissen von Null verschiedenen Residuen (16. 1) im Innern und von Residuen’ 
(16. 2) am Rande behaftet sein. Man kann aus deren Kleinheit nicht unmittel- 
bar auf die Giite der Approximation schliessen, wie manche Autoren anzuneh- 
men scheinen. Es kénnen sich vielmehr trotz winziger Residuen unter Um- 
standen Naherungslésungen einstellen, welche von der gesuchten Gitterlésung 
auffallig verschieden sind. Darauf hat schon F. WotrF [24] in ane 
Weise aufmerksam gemacht. | 

Zur Beurteilung der zu erwartenden Fehler gehen wir ahnlich wie S. GERSCH- 
GORIN [10] vor. Es mégen im Gittergebiet G’ u(x, y) die exakte und w(x, vy) 
die angenaherte Lésung des Gleichungssystems (4. 1) bezeichnen. Wir haben 
also unter Beniitzung der in (16. 1) und (16. 2) eingefiithrten Residuen : 


" 
q 


lu)=y, Ww)=yt+r; LU)=F, LU)=F+R. (229 


Daraus ergeben sich fiir die Differenz 


0'(x, y) = u'(x, vy) — u(x, y) (22. 2) 

die Relationen 
Ud'j=r, -L(6) = Rs (22. 3 
; 


Nun kénnen wir auf 6’ die in § 2 fiir andere Zwecke entwickelten Abschat- 
zungen anwenden. Setzen wir namlich 


r* Mam za meee a (22. 4) 
so haben wir an Stelle von (9. 8) und (12. 3) die Ungleichungen in 1(6’)| s7*9 
| L(6’)| S R*. Fithrt man ganz analog wie in (13.1) und (13. 2) die Hilfs- 


funktionen 2z(x, y) und w(x, y) ein, so gelangt man fiir eine hinreichend kleine 
Gitterseite 4 zu der (13. 3) entsprechenden Abschatzung 


| 6’(x, y) | S ma(x, y) . (22. 5) 
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Nur hat hier m an Stelle von (14. 2) der Ungleichung 


1 
Te fe. Max (2 7*, 2h R*) (225.6) 


zu gentigen, weshalb wir 


if 
Nii ore xy (2257) 


mit einer den maximalen Betragen der Residuen proportionalen, aber von h 
unabhangigen Konstanten 7, setzen diirfen. (22. 5) geht damit in 


F 1 
| o’(x, y) |S je 'M 2(x, ¥) (22. 8) 


uber, wo 2(x, vy) die Lésung der von A unabhangigen Randwertaufgabe (13. 1) 
bezeichnet. 

Man wird in einem Zahlenbeispiel die Ungleichung (22. 8) kaum fiir die 
Abschatzung von 0’(x, y) heranziehen. Einerseits wiirde sich die Miihe nicht 
lohnen. Anderseits ware 7 zu gross und daher die Schranke (22. 8) zu hoch. 
Denn beim Relaxationsverfahren kann man die Approximation trotz be- 
schrankter Stellenzahl erheblich weitertreiben, als dies die Abschatzung dieses 
_ Abschnitts erwarten liesse. Wie schon erwahnt, sorgt man ja dafiir, dass die 
verbleibenden minimalen Residuen von verschiedenen Vorzeichen sind, ja dass 
nirgends Blocke mit vorwiegend positiven oder negativen Residuen stehen 
bleiben. 

Dagegen kann man aus (22. 8) den Einfluss der Verfeinerung des Gitters 
auf den zu erwartenden Fehler 6’(x, y) erkennen. Er ist umgekehrt proportional 
zum Quadrat der Gitterseite. Soll also bei Verfeinerung des Gitters der durch Ab- 
_ rundung verursachte Fehler 0'(x, y) nicht vergrossert werden, so muss man gleich- 
zeitig die Anzahl der Dezimalen, mit denen man die Werte der Gitterfunktion 
schreibt, erhihen. Insbesondere ist beim Ubergang zu einem Gitter mit halber 
_ Gitterseite die relative Genauigkeit der Gitterwerte zu vervierfachen. 

Wenn man diese Vorsicht walten lasst, so k6nnen die Abrundungsfehler die 
Konvergenz bei Verfeinerung des Gitters nicht in Frage stellen. Deshalb diirfen 
wir das Differenzenverfahren mit Recht als abrundungsfest!) bezeichnen. 
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Summary 


Partial differential equations of the elliptic type are practically solved by 
finite-difference approximation. If the bowndary conditions involve derivatives, the 
approximation by finite differences can be a serious problem. In this paper it is 
shown how the normal derivative may be approximated for a curved boundary 
without any corner. The finite-difference solutions converge to the exact solution 
of the differential equation while the mesh-side tends to zero. This is proved 
under certain conditions. The method used is similar to that given by S. GERSCH- 
GORIN. The convergence rate is discussed in detail. The final part of the paper 
is concerned with the relaxation technique. It is proved that the residuals for inner 
and for boundary points converge to zero. Finally the rounding errors are discussed. 


{Eingegangen: 25. 10. 1951.) 


Zur numerischen Bestimmung periodischer Losungen 
- bei nichtlinearen Schwingungen 


Von LotHar CoLiatz, Hannover?) 


1. Die Randwertaufgabe 


Bei der Bestimmung periodischer Losungen bei nichtlinearen Schwingungen 
liegen nichtlineare Randwertaufgaben vor. Man hat fiir eine oder allgemeiner 
fiir ~ unbekannte Funktionen +,(¢),..., x,(¢) (etwa von einer unabhangigen 
Veradnderlichen ¢) » Differentialgleichungen, welche entweder ¢ explizit gar 
nicht enthalten mégen oder, falls ¢ explizit auftritt, in ¢ die Periode T besitzen 
mégen, und Randbedingungen der Periodizitat, etwa fiir Differentialgleichun- 
gen zweiter Ordnung, angeschrieben: 


X(t) = 4 (@+T), x(t) =2(¢+T). G=L..m) () 


Die Periode T des Vorganges wird haufig (zum Beispiel bei erzwungenen 
Schwingungen meistens) bekannt sein, andernfalls ist sie als weitere Unbekannte 
einzufiihren (vgl. Nr. 5). 

Will man mit einem Naherungsansatz eine angenadherte Lésung ermitteln, 
so erhalt man im allgemeinen fiir die beim Naherungsansatz noch verfiigbaren 
Parameter ein System von etwa q nichtlinearen Bestimmungsgleichungen. 
Dieses zu lésen, ist meist eine miihsame Aufgabe, besonders bei grdsseren 
Werten von g. Daher sind méglichst genaue Verfahren erwtinscht, bei denen 


1) Lehrstuhl fiir hhere Mathematik der Technischen Hochschule. 
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man schon mit wenigen Parametern und also mit einem nicht zu grossen nicht- 
linearen Gleichungssystem benutzbare Werte erhalt. Die folgenden Verfahren 
wird man jedoch nur anwenden, wenn die in der Differentialgleichung auf-— 
tretenden Funktionen und die Lésungsfunktion geniigend «glatt» sind und — 
stetige Ableitungen so hoher Ordnung wie die des unterdriickten Restgliedes 
besitzen. 


2. Kollokation mit Relaxation 


Zur Behandlung grésserer nichtlinearer Gleichungssysteme wird man oft” 
Iterationsverfahren benutzen; gute Erfolge zeigt dabei das Relaxationsverfah- 
ren. Durch einen vorlaufigen gréberen Naherungsansatz mit wenigen Para-— 
metern verschafft man sich einen Uberblick tiber die Lésungen und hat dann ~ 
bei der genaueren Rechnung bereits Ausgangswerte, die man iterativ ver- 
bessert. Das Verfahren seiim Zusammenhang mit der Kollokation, einem nahe- 
liegenden Naherungsverfahren, hier nochmals kurz beschrieben, da die Relaxa- 
tion auch noch im nachsten Abschnitt gebraucht wird. Es sei w,(¢) eine Nahe 
rungsfunktion fiir »,(¢). Fir w; macht man zum Beispiel einen Ansatz eines . 
endlichen trigonometrischen Polynoms: 


z ’ ’ 
x;(t) & w,(t) = DS (6; nose e+S, 52D sua t) (7 = 1, ..<, 9%)" ee | 


(mit so; = sp; = 0); fiir die 2p unbekannten Fourier-Koeffizienten c,, und 

s,; erhalt man 2 Bestimmungsgleichungen, indem man verlangt, dass die~ 
n Differentialgleichungen, wenn man w fiir x; einsetzt, an 2 Stellen, den- 
«Kollokationsstellen » 


10, hij eh ee 
2 p ’ 2 p ’ ’ 2 p ’ 
erfiillt sind. 
Beispiel: Es sei die Gleichung : 
x-+%+ (x +x’) =cost (3) 


(erzwungene Schwingung eines eingliedrigen gedampften Schwingers mit nicht- 
linearer Riickstellkraft) vorgelegt. Der Ansatz [unter Ausnutzung der Schief- 
symmetrie x(¢ + 2) = — x(t)] 

x(t) ~ w(t) = cy cost+ s,sin¢t (4) 
fiihrt mit den Kollokationsstellen 0 und a/2 auf 


Ss tce=1, —c¢,+s8=0 


7 0}=49) — a= . . . 
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mit der Losung 
Cp 0,960 5s, = 0)825 . (5) 


Bei dem Ansatz mit vier Parametern 
w(t) = cy cost + s,sint+ cs cos3¢+ s3sin3¢ (6) 


und den Kollokationsstellen t = 0, 2/4, 2/2, 3 2/4 erhalt man das Gleichungs- 
system 


S) Sts+ 3534 (cy + cs)8 Se abs 
Gt S,+ 56 Tsp Cie Sia Sa) == Sk 
7) 
—ty r St, - 8s, (Sy — S3)8 0), 0) 
(%—s,—11e, ise (—c, + 5, + ¢3 + $3)? = —1 
oder mit 
Cy + Cs =; €y + Sy = GE S35 = B; 
5, Ss =; G%+5,+¢,+ s,= 0; 
OL eae eae 
B2—2+ 44-864 12y—26=0, 
(8) 


y? = Cais B 44-3 0= 0% 
O-—-2—12Za+268+ 4y—8d=0. 


Das folgende Schema (Tabelle 1) zeigt die Durchfiihrung der Relaxations- 
rechnung. Es wird ausgegangen von den Werten 


Cp = UDO 8, 0S, i, = 0, sg 0S, 


die nach (5) durch die grébere Rechnung mit c, und s, nahegelegt werden; dann 
rechnet man jeweils die linken Seiten der Gleichungen (8), die «lunktions- 
werte», aus und stellt fest, wo sie noch am starksten vom Wert Null abweichen. 
Nun werden an den Ausgangswerten ¢,, $1, C3, $3 kleine Anderungen (im Schema 
eingeklammert angegeben) angebracht und die neuen Funktionswerte und die 
erfolgten «Anderungen der Funktionswerte» ausgerechnet. Diese Anderungen 
der Funktionswerte erlauben es, die Dosierung der Anderungen der Ausgangs- 
werte geeignet vorzunehmen, so dass die Funktionswerte sich dem Wert Null 
naihern. Ein grobes Gesamtmass fiir die Gitte der Naherung ist die im Schema 
als «Spanne» bezeichnete Summe der Betrage der Funktionswerte. 
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Tabelle 1 
o B y m) ot? Bp v3 3 A}) 
0,55 Lo 0,80 0;25)) 0,17 2,46 0,51 0,02 |—0,20 
(— 0,05) : 
0,55 oo 0,75 0525" | 0,27 2,46 0,42 0,02 |—0,15 
(- 0,07) 
0,48 dl 635) 0,75 0,25 One 2,46 0,42 0,02 0,13 
0,42 122 0,74 0,30 | 0,074 | 1,816 | 0,405 | 0,027 |—0,06 
0,408 | 1,206 0,741 | 0,314] 0,0679] 1,7540| 0,4069]| 0,0310|— 0,069 
(— 0,0001) (— 0,0001) : 
0,4079 | 1,206 0,7409 | 0,314] 0,0679] 1,7540| 0,4067| 0,0310|—0,0685 | —1,7536 
1) A=—1—4a4+3fB—y—0. *) B=—2+4a—8 B+ 12y—20. 
Ergebnis: 


Erste | Zweite | Dritte 
Naherung 
mit bis bis 
07,84 | base C5, Ss 


0,408 


3. Gewohnliches und verbessertes Differenzenverfahren 


Bei allen Arten von Differenzenverfahren legt man eine Intervallteilung 


ogee 


i A A ee a A MR Aly inlet 


a ee) 


zugrunde, etwa eine Teilung des Periodenintervalls in m gleiche Teile der 


Linge h = T/m und stellt Bestimmungsgleichungen fiir Gréssen x,, auf, welche 
Naherungswerte fiir x,(0 /) (fiir @= 0, 1, ..., m—1) sein sollen. Das gewéhn- 


liche Differenzenverfahren erscheint fiir den vorliegenden Zweck wenig geeig- 
net. Bei ihm werden die Differentialgleichungen fiir t = 9 h angeschrieben und 
dabei alle in den Differentialgleichungen vorkommenden Differentialquotienten 
durch Differenzenquotienten ersetzt. Bei grosser Maschenweite / und geringer 
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Tabelle 1 
C8) pt) Funktionswerte Anderungen E 
Ara? [2 pe C+? |D+ 63 der Funktionswerte eeaial 
F955 0,70 0,03 0.96 0,04 |—0,68 7A 
—0,35)_ | 0,90 0,02 |0,36 0,07 |-0,88 |0,05 |-0,60 0,11 |-0,20 11,33 


—0,42 -|- 0,06 0,24 |0,08 0 —0,04 |0,22 |-0,28 |-0,07 | 0,84 |0,36 
0,408 +-0,032 |-0,0011]0,0020-0,0011|-0,001 | | | _ {o,o052 
le ees 0,0006} 0,0004|— ees 0,0002) 0,0005/— 0,0016] 0,0001] 0,0008] 0,0022 


8) C=a+B—4y4+306. 4) D=2—1204+2B+4y—80. 


Anzahl von Unbekannten werden die erhaltenen Naherungen dabei oft viel zu 
grob (vergleiche das folgende Beispiel), und eine Maschenverfeinerung ist 
gerade bei den nichtlinearen Aufgaben wegen des Anwachsens der Unbekann- 
tenzahl ungiinstig. 

Als eine der verschiedenen mdglichen Verbesserungen des Differenzenver- 
fahrens ware das Mehrstellenverfahren') zu nennen. 

Wir beschranken uns bei der Beschreibung etwa auf eine unabhangige 
Funktion x(¢). Man verwendet Ausdriicke der Bauart 

i+p 
P= Ze, eI een one o (9) 

_ wobei bei einer fest gewadhlten k-ten Ableitung die a, und A, so bestimmt 
werden, dass bei Taylor-Entwicklung der Werte %;,) = «(vy h) und Ableitungen 
x{t) — x(y h) und damit des Ausdrucks P an der Stelle t;=7h die Ablei- 
tungen bis zu méglichst hoher Ordnung den Faktor Null erhalten, das heisst 
das Restglied von méglichst hoher Ordnung wird. Fiir die erste und zweite 
Ableitung zum Beispiel kann man so die Ausdriicke aufstellen 


%i-1) — X41 +z Been + 4 %H TF %41)) = Restglied 5.Ordnung _— (10) 


11) 


1) Vergleiche zum Beispiel die zusammenfassende Darstellung bei L. CoLLatz, Nuimertsche 
Behandlung von Differentialgleichungen (Springer, Berlin, Gottingen, Heidelberg 1951), S. 116. 
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Kommen in der Differentialgleichung zum Beispiel die Gréssen x, x, x vor, 


so hat man fiir jede der m Stellen t = @ h die Mehrstellengleichungen mit k = 1 ~ 
und 2 (wobei das Restglied durch Null und x, und xin durch Naherungen ~ 


%;, xa" zu ersetzen sind) und die Differentialgleichung anzuschreiben und er- 


halt so 3 m Gleichungen fiir die 3m Unbekannten ;, x;, x;. Im folgenden soll — 


eine Modifikation beschrieben werden, nach der man in vielen Fallen unmittel- 
bar auf ein Gleichungssystem von nur m Gleichungen fiir m Unbekannte ge- 
fiihrt wird. Diese Modifikation ist stets durchfiihrbar, wenn die Differential- 
gleichung sich auf die Gestalt bringen lasst: 


Lx} Sie a1 (12) 


wobei L[x] ein linearer Differentialausdruck in x und seinen Ableitungen mit 
konstanten Koeffizienten und N eine beliebige (nichtlineare) Funktion von ¢ 
und einer festen Ableitung x) ist. Viele bei nichtlinearen Schwingungsauf- 
gaben auftretende Differentialgleichungen lassen sich auf diese Form bringen. 
Dann verwendet man Ausdriicke der Gestalt 
itp 
P = pat (a, L[%,,)] + A v #3) (13) 
v=t—p 
und bestimmt wieder die a, und A, so, dass bei Taylor-Entwicklung von P 
an der Stelle ¢; = 7h das Restglied von méglichst hoher Ordnung wird. Bei 
Unterdriickung des Restgliedes und Einsetzen der Differentialgleichung erhalt 
man m Gleichungen fiir die # Unbekannten x), xl) ra wos: 
+p 
a: (a,N [x0 ¢]+ A, x0) =0. (¢=0,...,.m—1) (44) 
v=t—p 


Bei Schwingungsaufgaben werden fiir L[x] oft Linearkombinationen aus 


ean ae 


x und x oder x und x in Betracht kommen; daher seien hier als vielleicht wich- ~ 


tigste Fille die beiden Formeln (15) und (18) angegeben, wobei wir 7 = 0 wah- 

ie Ausrechnung der Taylor-Entwicklungen kann hier unterdriickt wer- 
Fir I[x]= x + s x (beliebige Konstante s) gilt mit sh =o 

(72 — 360+ 3.0%) x_, — 144 %)+ (72+ 360 — 30°) x, | 

+ h®{(— 6+ 30+ 0%) L[x]_, + (— 60+ 40%) L[x], (15) 

+ (— 6— 36 + 6%) L[x],}= Restglied 5. Ordnung. | 


Diese Formel ist geeignet fiir Schwinger mit linearer Dampfung und _nicht- 
linearer Riickstellkraft, also fiir Differentialgleichungen 


x-++s%*% + f(x) = Pit). 
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Pur s > oo und s = 0 sind hier die obigen Formeln (10) und (11) als Spezial- 
falle enthalten [allerdings zeigt sich dabei in (11), dass in der Taylor-Entwick- 
lung die Glieder 5. Ordnung fortfallen und das Restglied von 6. Ordnung 
wird]. Fiir o = 2 ergibt sich die nachher beim Beispiel gebrauchte Formel 


© (%_1 — 6 % +5 x) +h? {L[x]_,— 11 L[x],—2 [x], }= Restglied 5. Ordnung. 


(16) 
Fir L[x] = x + s x lautet eine genauere Formel: 
(3 o* — Il o? — 45 6 + 45) x_, + (84 0? — 16 o? — 3600 + 720) x_, 
+ (5407 — 1530) x) + (—84 o? — 16 62+ 3600+ 720) x, | 

(17) 


pet o — 116" + 450+ 45) x, -+ 12h? { (3 o? + 50 — 10) Lx], 


+ (96% — 55) L[ x] + Bo? — 50 — 10) L[ x], }= Restglied 7. Ordnung 
mit den Spezialfallen 
o=0: 3 (%.+16%_, — 34%) +16 %,4+ x) —4h? (2%_,+11% 4+ 2 x,) 

= Restglied 8. Ordnung , 

S>0o: % 5+ 28 %_,— 28 x — % + 12h (x. +3 x%9)+ %4) 
= Restglied 7. Ordnung. 

_ Ferner sei genannt fiir L[x] = * + s x 

(—12+4sh?) x) + (6+ sh?) (4_14+ x) | ie 

+ 3hf («+ s x)_,—(¥ +s x),}= Restglied 5. Ordnung. | 


Diese Formel ist brauchbar bei Schwingern mit nichtlinearer Dampfung und 
linearer Riickstellkraft : 

x + g(x) +sx= Pit). 
Fiir s > oo erhalt man wieder die obige Formel (10), fiir s = 0 folgt 


2 (%_4 — 2 %9 + %) + h (%_ — %) = Restglied 5. Ordnung. 
Die vieluntersuchte Gleichung 
x + f(x) %+c¢% = rt) (19) 


hat nicht unmittelbar die Gestalt (12), lasst sich aber auf diese Gestalt um- 


schreiben, indem man ; 


x(t) =4(t) und x j(x) dx = F(x), [re di = Rit) 


x ty 
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setzt. Die Gleichung 

¢+ 5 (F@)+c#=r00) (20) 
lasst sich einmal integrieren: 

Z+F()+ez2=C+ Ri), (21) 


wobei die Konstante C belanglos ist, da eine additive Konstante bei z die 
Funktion x(¢) nicht beeinflusst; man kann also zum Beispiel C = 0 setzen. 
Gleichung (21) hat nun die Form (12) und kann mit der Formel (18) behandelt 
werden; in dieser Formel treten dann die Gréssen 2, also die interessierenden 
Werte x, aber nicht z selbst auf. 


4. Beispiel 


Es werde die Anwendung auf ein System beschrieben: 
x+2x4+f(*)+x-—y=0, 2j+y—x=cosst (22) 


(Schwingungsanordnung nach Figur 1). Die Riickstellkraft werde hier durch 


Fig. 1 


Schematische Schwingungsanordnung zu dem Beispiel in Nr. 4. 


{(x) = x + @ x* festgelegt, wobei fiir @ = 0 eine lineare Aufgabe entsteht mit 
der periodischen Lésung: 


x= 0,0055 cos t — 0,0805 sin = 2, | 


(23) 
y = — 0,2555 cos t+ 0,0205 sin t. | 


Wir legen der Rechnung die (noch sehr grobe) Maschenweite h = 1 zugrunde und 
haben dann fiir x und y jeweils vier unbekannte Funktionswerte, die wir hier 
kurz mit a, b,..., 6 bezeichnen: 


ie oa” ae see al 
Na OF RCM INLETS 
J: 8 ya Ome 


me 
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Das gewohnliche Differenzenverfahren fiihrt dann auf acht Gleichungen, wobei 
wir von jedem Satz von vier Gleichungen nur die zu ¢ = 0 gehorige Gleichung 
anschreiben: 


an ah ay (ae 0 a2) og 0 = 0, 
(-404+284+26)+a-—a = 1. 


Aus den letzten vier linearen Gleichungen driicken wir «, B, y, 6 durch a, b, c, d 
aus: 


2Zla=—7+a+664+ 8c+ 6d, 


218 = Ga bE 6c 8 a. 
2ly=7+8a+66b+ c+6d, 
21 6 = Oa Sa Oede @ 


und erhalten entsprechend vier Gleichungen fiir a, b, c, d: 


2loea*+7—. a+36b 8c Girne 100, 
21 010" 64a b+ 36¢ 8d=0, 
2Zloc?—7 8a 66 e=— F560 =O. 
21 0 a? +36a— 8b— 6c— d=0. 


Da wir die Ergebnisse mit denen des genaueren Verfahrens vergleichen wollen, 
fiihren wir gleich noch die Rechnung mit Benutzung der Formeln (15) und (18) 
vor. Auch hier hat man zwei Satze von je vier Gleichungen, von denen wir 
jeweils die erste Gleichung hinschreiben: . 


— 364+ 3064+ 6d—1la+ 2244+ llga®—2f+40d 
+200?+6—2d—oda*=0, 
oder 0 (11 a2 + 262— d?) —14a+ 346+4d—11a—2f8+6=0, 
und —48a+ 246+ 246—10a+10a—10—6+f8—d+6=0. 


Wieder driickt man zuniichst die «, B, y, 6 durch a, 0, c, d aus: 


209 a = —55+ 20a+ 5764+ 75c+ 574, 


209 B = 57a+ 206+ 57c+ 75d, (24) 
209y= 554+ 7544 576+ 20¢+ 5/4, 
AC oe 57a+.756+ 57c+ 20d 


und erhalt einsetzend vier Gleichungen fiir a, D, c, d: 
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209 9 (11 a? + 2.52 — d?) + 605 — 3203 a+ 6514 b — 882c4— 77Dad—= 08 
209 0 ( —a?+110?+4+ 2c?)— 165+ 79 a — 3203 6-= 65146 = -c820— ne 
209 0 ( — 624+ 11¢c?2+ 2d?) — 605 — 882a+- 796—3203 c+ 65144 ¢d=—08 
209 0 ( 2a? — c2+ 11d?) + 165 + 6514 a — 8826+ 79\¢— 3203 cd=07 


es empfiehlt sich fiir die folgende Relaxation, diese Gleichungen nach a*, bes 
c2, d* aufzuldsen. Dann erhalt man 


13585 9 a+ 3575 — 15512 a+ 40185 6 — 12399c+ 1311d=0, 


13585 0 b? 4 1311a@— 15512 6+ 40185 c — 12399d=0, (25) 
13585 9 c? — 3575 — 12399a+ 1311 6— 15512c+ 40185d=—0, 
13585 @ d? + 40185 a— 12399b+ 1311¢—15512d=0. 


Es wurde zunachst der lineare Fall (90 = 0) durchgerechnet und die Werte mit 
der exakten Lésung verglichen. Die Tabelle 


Linearer Fall (9 = 0) 
(Absolute Fehler in Klammern beigefiigt) 


Gew@éhnliches Verbessertes Differenzenverfahren 
Exakter Wert 


Differenzenverfahren nach (16) nach (17) 


0,00547 | —0,0270 (— ( 0,0018) | 0,0073 ( 0,0018) 
—0,08053 | —0,1622 (~0,0817) | —0,0914 (—0,0109) | —0,0824 (—0,0019) 
—0,25553 | —0,3243 (-0,0688) | —0,2651 (—0,0096) | —0,2577 (—0,0022) 

0,02047 | 0,0541 ( 0,0336) | 0,0240 ( 0,0035) | 0,0211 ( 0,0007) 


0,0324) 0,0073 


Tabelle 2 


0 —0,1 0 0,1 0 135,85 0 135,85 |—312,4 3Tt3 
0 —0,1 —0,01 0,1 0 135,85 1,4 |135,85 |—188,4 — 90,6 
0 —0,1 —0,01 0,09 0 135,85 1,4 |110,0 |—201,5 33,4 
0,002 |—0,095}]—0,013} 0,088 0,05 }122,5 2,59 | L0S50 3,0 |—137,3 


0,953 |— 1,055 | —1,042 |— 0,940} 12338 |15120 |14750 |12004 |-12350 |-15102 


(Die zweite Lésung) 


1) A = 8575 — 15512a +---. %) B= 13lla+---, 


EEE 


<eeeeel te. 6 TTT 


2 a er Me 


i A eet a, MER A ge me i At 


Vol. III, 1952 Zur numerischen Bestimmung periodischer Lésungen 


203 


zeigt, dass dabei das gewohnliche Differenzenverfahren recht schlechte Werte 
liefert und dass die Mehrarbeit bei Aufstellung der Gleichungen (16) sich gelohnt 
hat. Man wird daher wohl auch im nichtlinearen Fall zu den Formeln (15) und 
(18) mehr Vertrauen haben diirfen als zum gewohnlichen Differenzenverfahren. 
Es ergeben sich zwei periodische Lésungen, das Schema zeigt fiir 9 = 1 die 
Durchfiihrung der Relaxation fiir die beiden Lésungen. Die Relaxation arbeitet 
sehr gut, wenn man sich erst einmal in das Verfahren hineingedacht hat. Da 
die Anlage der Rechnung ganz entsprechend wie beim vorigen Beispiel ist, 
ertibrigt sich eine genauere Beschreibung (siehe Tabelle 2). 
Ergebnis der Rechnung: 


Erste Lésung Zweite Losung 


ai mes? digo lst. 9 
0 0,002 — 0,269 = 0,953 —1,272 
1 — 0,095 0,020 —1,055 — 0,982 
2 — 0,013 0,260 =1,042 —0,721 
3 0,088 =0,029 —1,012 


5. Beispiel: Periode unbekannt 


Ist die Periode T nicht bekannt, so hat man sie als weitere Unbekannte ein- 
zufiihren. Dafiir kann man, wenn in den Differentialgleichungen die Zeit ¢ nicht 
explizit auftritt, den Zeitpunkt ¢ = 0 beliebig festlegen, zum Beispiel in einem 
Augenblick maximaler Auslenkung einer Grisse etwa %,(0) = 0 wahlen. Ledig- 
lich zur Erlauterung der Methode werde folgendes einfache Beispiel behandelt: 


£4 B%—#3)+%=0, 


Tabelle 2 
Funktionswerte Anderungen eee 
©) D") 5 5 der Funktionswerte WY 
A+ NaIB+Nb?2\C+ Nc7|D+ Nd? 

312.4 |-311,;3 | —312 447 only —175 1246 
467,5 |-324,4 | —188,4 45,3 468,9 |—188,5| 124 |-402 Pap Nes) 890 
65,6 }-169,3 | —201,5 | 169,3 67 — 59,3 | —13 124 |-402 WS 495 
13,6 |}-123,8 3,0 } -14,8 15,9 |— 18,8 53 
-14752|+-12000| —11 18 —3 4+ 36 


8) C =—3575— 12399a + -«-. 


4) D = 40185a4+---. 


. 
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welches so einfach ist, dass man es fiir technische Zwecke sonst wohl durch 
Zuriickfiihrung auf eine Differentialgleichung 1. Ordnung und graphisch mit 
Hilfe der Liénardschen Konstruktion behandeln wiirde. Es werde also 
x(0) = x(T) = 0 gesetzt. Bei der ganz groben Maschenweite 4 = 7/4 hat man 
fiir die beiden Unbekannten h und x(k) = a nach (18) die beiden Gleichungen 


6h(—3a+a3)=0, (—12+4h*2)a=0 


mit der Lésung (die triviale Lésung # = a= 0 lassen wir ausser Betracht) * 


h=a=V3 und T=4h x 6,93. Diese Werte werden als Ausgangswerte fiir 
die Relaxation bei der kleineren Maschenweite h = T/12 benutzt. Fiir die sechs 
Unbekannten h, x(h) = a, x(2h) = 6, x(3h)=c, x(4h) =d, x(5h) =e (vel. 
Figur 2) lauten die Gleichungen nach (18) 


PNPTA ay J 42T § 


Fig. 2 
Verschiedene Naherungen zu dem Beispiel in Nr. 5. 


3h[—3e+ e—3a+a'| + (6+ h?)(—e+a)=0, 
3h —3b+ b5]+ (—124 4h%) a+ (64h) ( b) =0, 
3h[ 3a—a®—3c + c8]+ (—124+ 4h) b+ (64+4%)( atc)=0, 
3h[ 3b—b3—-3d+4 d*]+ (—-124+4h%)c+ (6+h2)( 64a) =0, 
3h[ 3c—c8—3e+ e]4+ (—12+4h%) d+ (6+h2)( c+e)=0, 
3h[ 3d—d3 ]4(C124+4A ec HOT AM 2 Oe Oe 


Im Relaxationsschema wird man h, b, d als drei Unbekannte einfiihren. Aus 
ihnen berechnen sich aus drei der angegebenen sechs Gleichungen die Gréssen 
a, c, e. Zu den jeweils gewihlten Werten von h, b, d bestimmt man die zugehé- 
rigen Werte a, c, e und stellt fest, wie weit die linken Seiten der restlichen drei 
Gleichungen dann die Zahl Null annahern. Durch Anbringen kleiner Variatio- 
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nen bei h, b, d gelangt man schnell zu einer hinreichend genauen Naherungs- 
lésung: 


Pa U,000, a= 1/5, U=202, €= 1,63, d= 1,54, ¢=1,19. 


Figur 2 zeigt den Verlauf dieser Naherungslésung und zum Vergleich gestrichelt 
auch die grobe Naherung mit h = 7/4. 


Summary 


In order to solve a periodic problem, the solution of its boundary value 
problem is needed. We ask for a method with a high order of accuracy and with 
nonlinear equations containing only a few unknowns. In this case the use of 
ordinary finite differences is not recommended. The ‘‘Mehrstellenverfahren”’ uses 
terms (9) and is always applicable, but gives a nonlinear system with a large 
number of unknowns. The method described in this paper is applicable for 
equations (12) and uses the terms (13). A, and a, are evaluated from the fact 
that the remaining terms of Taylor’s expansion are of an order as high as 
possible. Equations (15), (16), (17), and (18) are examples of such formulas. Some 
numerical examples are also given. 


(Eingegangen: 29. 11. 1951.) 


Schallgeschwindigkeit in pordsen Medien 


Von HEINRICH LABHART, Ziirich?) 


1. Einleitung 


Die Kenntnis der Abhangigkeit der Schallgeschwindigkeit von der Porositat 
fester K6érper ist von vielseitigem Interesse. Einerseits ist es dadurch méglich, 
durch rasche, serienmaBig durchfiihrbare Messung der Schallgeschwindigkeit 
auf das Porenvolumen von Festkérpern, insbesondere von Gesteinen und kera- 
mischen Massen verschiedenster Art, zu schlieBen. 

Andrerseits ist die Kenntnis der Schallgeschwindigkeit unerlaBlich fiir die 
Dimensionierung der heute immer mehr an Stelle von Einkristallen als elektro- 
mechanische oder magnetomechanische Schwinger verwendeten keramischen 
Sinterkérper, wie BaTiO, oder magnetostriktive Ferrite. 

Kiirzlich hat GASSMANN®) den EinfluB einer Porenfiillung durch eine kom- 
pressible Fliissigkeit auf die elastischen Eigenschaften eines pordsen Mediums 
behandelt. Uber den Zusammenhang zwischen dem Elastizitatsmodul des po- 
résen Materials und demjenigen der homogenen Geriistsubstanz liegen bisher 


1) Abteilung fiir industrielle Forschung des Institutes fiir technische Physik an der ETH. 
2) F.GassmAnn, Vjschr. naturf. Ges. Ziirich 96, Heft 1 (1951). 
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nur die Untersuchungen von SgurreE!) und BARTELS?) an Sintereisen, von 
RyYSCHKEWITSCH) an keramischen Massen und die erst wahrend der Fertig- 
stellung der vorliegenden Arbeit erschienenen Versuche an Ni-Zn-Ferrit von 
WEIL) vor. In diesen Arbeiten finden sich jedoch keinerlei Ansatze fiir eine 
theoretische Deutung des Effektes. 


2. Experimentelle Untersuchungen und Ergebnisse 


Fir die vorliegenden Experimente wurden magnetostriktive Ferrite der 
angendherten stéchiometrischen Zusammensetzung Nig,, ZNo,> Fe,O, verwendet. 
Dieses Material laBt sich durch gemeinsames thermisches Zersetzen der zusam- 
men gemahlenen Metalloxalate als feines Pulver erhalten, das mit einem orga- 
nischen Binder versetzt und in einer Matrize zu Stabchen verpreBt wird. Durch 
Sintern dieser Stabchen bei Temperaturen zwischen 1100°C und 1280°C erhalt 
man Porenvolumina von 26% bis 4°. Die fertigen Stabchen haben eine Linge 
von etwa 8cm und einen anndhernd quadratischen Querschnitt von etwa 
0,3. cm?. Fiir die Messung der Resonanzfrequenz werden die Stabe in einer 
zylindrischen Spule, genau in der Mitte, wo bei Anregung in der Grundfrequenz 
ein Schwingungsknoten liegt, mit vier spitzen Schrauben festgeklemmt. An 
allen anderen Stellen schwingen sie frei. Dem erregenden Wechselstrom in der 
Spule wird zur Vormagnetisierung ein Gleichstrom iiberlagert. Da die Schwin- 
gungsgiite Q dieser Stabe fast immer erheblich gréBer als 1000 ist, zeichnet sich 
die Resonanzfrequenz durch ein sehr scharfes Maximum der Impedanz aus 
und kann genau festgelegt werden. Schwingt ein Stab der Lange / in seiner 
tiefsten Eigenfrequenz, so berechnet sich die Schallgeschwindigkeit zu 


[gee ee 


Die Eigenfrequenz der verwendeten Stabe liegt zwischen 25 und 45 kc. Die 
durch die Vormagnetisierung bedingte Anderung des Elastizitaétsmoduls ist 
groBenordnungsmabig kleiner als der EinfluB des Porenvolumens und muB 
daher nicht beriicksichtigt werden. Die Bestimmung des Porenvolumens ge- 
schieht pyknometrisch. 

In Figur 1 ist die Schallgeschwindigkeit von 30 Stabchen von annahernd 
gleicher chemischer Zusammensetzung in Funktion des Porenvolumens dar- 
gestellt. 

Mit einer gewissen durch Inhomogenitaten und kleine Fehler der Stabchen 
bedingten Streuung scheint zwischen Porenvolumen und Schallgeschwindigkeit 
eine eindeutige Abhangigkeit zu bestehen. Aus der Schallgeschwindigkeit 


1) A. Squire, Powder Metallurgy (Mapleton House, New York, 1947) 
2) H.J. BARTELS, Mitt. Kongr. Pulvermetallurgie Graz 1948. 
3) A.RyscuKewitscu, Ber. dtsch. keram. Ges. 25, 95 (1944) 
4) L.Weit, Ann. Inst. Fourier 2, 207 (1951). 
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schlieBt man fiir einen Ferrit mit einem Porenvolumen von 10° und einer 
Dichte von 4,80 g/cm® auf einen Elastizitatsmodul EF = 1,46 - 1012 Dyn/cm2, 
der mit demjenigen von Stahl (E = 2,2 - 1012 Dyn/cm2) vergleichbar ist. 


3. Theoretische Deutung 


Der Spannungszustand eines pordésen K6rpers ist auBerordentlich kompli- 
ziert. Es ist daher unmoglich, durch eine strenge Rechnung aus den Elastizi- 
tatskonstanten des Gertistmaterials diejenigen des pordésen Stoffes zu erhalten. 

Deshalb soll eine Naherung versucht werden, welche in Folgendem besteht: 

Man betrachte einen pordsen Zylinder von 1 em? Querschnitt, der unter 
dem Zug Z steht (Figur 2). Die Poren sind durch Kugeln schematisiert. Legt 
man senkrecht zur Zugrichtung eine Schnittebene, so wird diese einige Poren 
schneiden und nahe bei anderen vorbeigehen. Die Summe der Projektionen der 
geschnittenen und der Ebene unmittelbar benachbarten Poren auf die Schnitt- 
ebene werde als Porenflache ~, bezeichnet. Der Zug Z verteilt sich somit auf 


Zé 


v-10-5 cm sekt 


Schniltebene 


Fig. 1 
Schallgeschwindigkeit in Abhangigkeit vom Porenvolumen. 


e MeBpunkte; — theoretische Kurven. 


die Flache (1 — #,) und bewirkt bei einem Elastizitatsmodul Fy des Geriist- 
materials eine relative Dehnung des Zylinders um 


yA Ne ie Se 
ae) Be 
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Der effektive Elastizititsmodul FE des porésen Kérpers betragt demnach 
E=E,(1— p;). 


Sind 0 die mittlere Dichte des pordsen Kérpers, gg die Dichte der Geriistsub- 
stanz und #, das relative Porenvolumen, also 


0=0o(1— fy»), 


so gilt fiir die Schallgeschwindigkeit in geniigend schlanken Staben 


v= |= - 2.428 = 4% 
0 Qo 1-— pf, Sat — Pe 
Wo vy die Schallgeschwindigkeit in der homogenen Geriistsubstanz bedeutet. 

Es besteht nun die Aufgabe, die Porenflache #; aus dem Porenvolumen , 
zu berechnen. Offenbar darf man #; nicht einfach gleich dem Mittelwert der 
Schnittflache der betrachteten Ebene mit den Poren setzen. Auf diese Weise 
wiirde man ~; = ~, und somit eine vom Porenvolumen unabhangige Schall- 
geschwindigkeit v = v, erhalten, was durch das Experiment absolut nicht be- 
statigt ist. Diese Diskrepanz riihrt natiirlich davon her, daB dem komplizierten 
Spannungszustand, der sich in der Umgebung einer Pore ausbildet, auf keiner- 
lei Weise Rechnung getragen wurde. Mit anderen Worten: Die Wirkung einer 
Pore erstreckt sich nicht nur auf ihr Volumen, sondern auch auf ihre nahere 
Umgebung. 

Im folgenden wird versucht, diesem Umstand Rechnung zu tragen. 

Dazu soll angenommen werden, daB die Poren einzelne, im allgemeinen 
gegeneinander abgeschlossene Hohlraume bilden, deren Gestalt nicht allzu- 
stark von der Kugelform abweiche. Ferner soll jede Pore von einer fiktiven 
Wirkungskugel umgeben gedacht sein (vergleiche Figur 2). Der Radius R; 
dieser Wirkungskugel sei so gewahlt, daB das Volumen der Wirkungskugel im 
selben Verhaltnis zum Volumen der Pore steht wie das Gesamtvolumen des 
K6rpers zum Gesamtporenvolumen. Liegt der Mittelpunkt einer Pore naher 
als Rk; bei der Schnittebene, so ist die Projektion dieser Pore als Beitrag zur 
Porenflache mitzuzahlen. Alle Poren, die weiter weg von der Schnittebene 
liegen, werden vernachlassigt. 

Die Poren seien entsprechend ihrem Volumen in verschiedene Sorten 
(Index 7) eingeteilt. Existieren von einer Sorte x; Poren pro Kubikzentimeter 
und ist der Beitrag dieser Sorte zum Gesamtporenvolumen #,,, so betragt 
offenbar das Volumen einer Pore dieser Sorte #,,/1;. Die Projektion o; dieser 
kugelig angenommenen Pore auf die Schnittflache wird somit 


of 3\28 [ Pip \2? 
o,= 78 (=) (es) . (1) 
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Der Radius R, der Wirkungskugel berechnet sich unter Beriicksichtigung 
der oben formulierten Bedingungen 


4a Rr? 4a R3 a Pin 1 


yay a Se ———— 2 
ay 5 Ha | und 


Fi 3 Z P. v 


3} Pi, 1/3 
is (2, ee) 


Die gesuchte Porenflaiche , wird somit 


3 P;, 1/3 , 3 \2/3 P; 2/3 
id et ape PAs : High jose eee 
ge Rae: Sas] nz? (Z) il 


a 


Zu 


3 2’ Pig 3 2/3 


Be Db ane 


Die Schallgeschwindigkeit wird damit 
‘1 — (3/2) p28 
v= |/ vee bis 


Diese Ableitung kann leicht in zweierlei Beziehung verallgemeinert werden. 

1. Fiir nicht kugelige Poren ist die mittlere Projektion o; der Poren der 
Sorte 7 nicht mehr nach Gleichung (1) zu berechnen, sondern sie wird um 
einen von der Porenform abhangigen Faktor 6 groBer. 


- 3 \2/3 ( Pin \?/8 
a 1/3 —_— 2 
Me y ie ( a ( nN; : 


Fiir ein Parallelepiped der Kantenlangen a, b, c berechnet sich 6 beispiels- 
weise Zu 


5 if (3 Ie (ab + be + @ Cc) 
eee NS 2 (abc)? 
Man erhalt somit 
Pardee sr et ones «ee emg On 00: 
fiir eitien!. Wirlelateecs ae 4. Ges eae) Omer LiA4: 
fir @ 0 Sl ea eS ee a Oe 16; 
CT ee a 20 ee ei a 0 2, oO; 


2. Durch die Zugrichtung ist im Material eine Richtung ausgezeichnet. Der 
einzufiihrende Wirkungsbereich einer Pore ist daher nicht notwendig kugel- 
symmetrisch, sondern méglicherweise in Zugrichtung gedehnt oder zusammen- 
gedriickt. Um dieser Méglichkeit Rechnung zu tragen, kénnen Rotationsellip- 
soide als Wirkungsbereiche eingefiihrt werden. Dabei soll immer angenommen 


ZAMP III/14 
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werden, daB die Lineardimension jeder Pore klein gegen die Lineardimension 
des zugehérigen Wirkungsbereiches sei und daB deshalb die Form der Pore die 
Form des Wirkungsbereiches nicht beeinflusse. Unter den gleichen Nebenbe- 
dingungen wie vorher erhalt man fiir die in Zugrichtung liegende Halbachse 0; 


wo y gleich dem Verhaltnis der Halbachse in Zugrichtung zu den zwei gleichen 
Halbachsen senkrecht zur Zugrichtung bedeutet. Fiir y > 1 ist also das Ellip- 
soid gestreckt, fiir y < 1 dagegen abgeplattet. 

Fiihrt man die Berechnung der Schallgeschwindigkeit wie vorher, jedoch 
unter Beniitzung von o; und 4; statt o, und R; durch, so ergibt sich 


1 — (3/2) 6 y2!3 p23 
v= %|/ oe 3A 


: 3 2/3 
oder mit m = — by"! 


v=%, +e 


4. Vergleich der Experimente mit der theoretischen Formel 


In Figur 1 sind unter Annahme von Vy = 6,18 - 10° cm/s die theoretischen 
Kurven fiir m = 1,5, m = 1,4 und m = 1,3 mit den MeBpunkten von 30 Staben 
zusammen eingetragen. 

Man sieht, daB diese einfache Theorie die GréBe des beobachteten Effektes 
durchaus richtig trifft. Fiir kleine Porenvolumen liegen jedoch die Werte naher 
an der Kurve fiir m = 1,3, wahrend fiir groBe Porenvolumen die Ubereinstim- 
mung mit der Kurve fiir m = 1,5 besser ist. Dies kann dadurch erklart werden, 
daB bei kleinen Porenvolumen die Sintertemperatur héher war und somit die 
Poren im allgemeinen kugeliger wurden. Dadurch wird 6 und somit auch m 
kleiner als bei den groBen Porenvolumen, wo infolge der niedrigen Sintertem- 
peratur die Poren starker von der Kugelgestalt abweichen. 

Es ist bemerkenswert, daB die Werte von 6 und y beide nahe bei eins liegen 
und somit eher als Korrekturen wirken. 

Setzt man namlich unter Annahme von Kugelporen bei kleinen Poren- 
volumen 6 = 1, so erhalt man sofort aus m = 1,3 fiir das Achsenverhiltnis der 
etwas abgeplatteten Wirkungsellipsoide y = 0,81. Mit diesem Wert fiir y erhalt 
man bei groBen Porenvolumen 6 = 1,15, was Poren entsprechen diirfte, deren 
Gestalt wenig von der Kugelform abweicht. 

Der Wert fiir den Elastizitatsmodul Ey des porenfreien Ferrits ergibt 
sich aus unserer theoretischen Extrapolation mit V, = 6,18 - 105 cm/s bei An- 


a 


Ag ne 
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nahme der aus den Réntgendaten berechneten Dichte 09 = 5:33 g/cm? zu 
2g = 2,04 - 1012 Dyn/cm?, was mit dem Wert 2,0 - 1012 Dyn/cm2 von WEIL in 
ausgezeichneter Ubereinstimmung steht. 


5. Diskussion 


Aus diesen Ausfitthrungen geht hervor, daB sich unter den gemachten, sehr 
rohen Annahmen (Querschnittsverringerung durch Porenflache und Einfiih- 
rung von Wirkungsbereichen der einzelnen Poren) die experimentell gefundene 
Abhangigkeit der Schallgeschwindigkeit vom Porenvolumen recht gut theo- 
retisch wiedergeben laBt. Es ist aber nach den bisherigen Uberlegungen auch 
zu erwarten, daB die Schallgeschwindigkeit nicht nur vom Porenvolumen, son- 
-dern auch, allerdings in geringerem AusmaB, von der Porenform und wahr- 
scheinlich von der GréBenverteilung der Poren abhangt. 

Daher ware es interessant, Versuchsreihen mit méglichst konstanter Poren- 
form und variablem Porenvolumen in der Weise durchzufihren, daB alle Stabe 
bei der gleichen Temperatur gesintert, aber mit verschiedenem PreBdruck ge- 
formt werden. — 

Weitere Daten iiber das elastische Verhalten kénnten aus der Messung von 
Biegeschwingungen an diesen Staben gewonnen werden. 

Durch oszillographische Messung der Abklingdauer frei schwingender Fer- 
ritstabe haben wir die innere Dampfung ermittelt. Unsere Versuche lassen noch 
keine einheitliche Interpretation zu. Die Verhaltnisse werden durch das Auf- 
treten von magnetomechanischer Hysterese kompliziert. 


Summary 


Measurements made show the very close relation between sound velocity and 
pore volume in a Ni-Zn ferrite of fixed composition. Schematizing the compli- 
cated state of stress and strain in a porous medium ina rough manner, a theoretical 
expression for the sound velocity is derived that fits the experimental curve 


astonishingly well. 


(Eingegangen: 25. 10. 1951.) 
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Bending of a Thin Elliptic Plate of an Orthotropic Material 


under Uniform Lateral Load 
By Yosro Onasi, Tokyo, Japan?) 

; 
In this paper, the bending problem of a thin elliptic plate of an aeolotropic — 
material subjected to a uniform load has been treated. The material of which 
is assumed to have two orthogonal axes of symmetry with respect to its elastic — 
properties. Recently the solution for a circular plate was obtained by OKUBO?). — 
In the case treated by him, however, the direction of the symmetrical axes of © 
the material are out of the question. While in cases of an elliptic plate, the © 
directions of the symmetrical axes come into question, because of the unsym- — 
metrical shape of the plate with respect to its center. 

We shall assume the plate to be of uniform thickness /, and shall take 
x-, y-axes in the plane of the plate before bending, the directions of which are — 
chosen in the directions of the symmetrical axes of the material. Then the — 
deflection of the plate satisfies the differential equation’) . 


04w 04w 04w 


dD, Ox! + 2 (D, + Dy) Ox? ov? 3 Oye q; (1) 


where gq is the intensity of a uniformly distributed load, S,, Sy., Soo, and Sg. 
are elastic constants defined by 


Caw = S11 Oy + Sie Oy, Cyy = Sy9'Fx + Sop Gy, Coy = Seg Torys 


and 
Ae Ss h?S : 
D: x Sis Be Se SBE L Lara 
a 12 (Sy, Sop Sig) id 12 (Sq, Seo — Sie) ‘ 
Aes hs 
Dy = fine Wiese. 


If x’, y’, and # represent the axes of the ellipse and the angle between 
x and x’, respectively, then 


x’ +49’ =e (x 4+ ty). (2) 
1) The Railway Technical Laboratory. 
*) H. Oxuso, J. appl. Phys. 20, 1151 (1949). 
3) S. TimosHENKO, Theory of Plates and Shells (McGraw-Hill Book Company, Inc., New York, 
1940), p. 189. 
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Transform the Cartesian coordinates x’, y’ into the curvilinear coordinates by 
Ro 19 == cosh(« + 7B), (3) 


then « = const represents a family of confocal ellipses and B = const represents 
a family of hyperbolas orthogonal to the former. Assuming that « = a» rep- 
resents the periphery of the plate, we have 


i COsmic, Ol C Sink io), (4) 


where a and 0 represent the half lengths of the major and minor axes of the 
elliptic plate, respectively. 

The general solution of the differential equation (1) can be written in the 
form!) 


W= Wt Rif(x+1h, y) + fo(x +7 ke y)I, (5) 


where W, is a particular solution of (1), /,, /, are arbitrary functions, R denote 
the real part of a complex function and 


ky = == = = = pe et ——= = saa ART , 
VD, =F D, =i {(De a D,)? _ D, Dee 


(a == = == ee 


1. Built-in Edge 


First, we shall consider an elliptic plate clamped at its periphery. We shall 
take for deflection the expression 


rs We \ 
w=4 (47+ ie (6) 


a? ge 


_ We see that w satisfies the boundary conditions w = 0w/d« = 0, and also the 
differential equation (1) provided 


Cl 
His L, cos*# + 2 L, cos?6 sin?é + Lz sint® ’ 
2D MIG Da. 24D 
L, Sea a2 b2 ate bt? 


(7) 
24 D 3 
Ly = 2458 + 8 (D+ Dy (Sr -— Gar t a) + aoe 


24 D, 16 (D, + D,) 24 Ds 
L3;= b4 a2 b2 a4 


1) H. Oxvuzo, J. appl. Phys. 20, 1151 (1949). 
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If be the angle which the normal to the curve « = const makes with the 
axes of x’, the normal makes an angle y — # with the axis of x, and so the 
bending moments can be written in the forms 


1 0?w 
EO =o Dyce Doi a) cos 2 (p — 3) Ox? 
1 7 0?w 
= 3 {Dats (Dz — Ds) cos2 (p — 8) Ov? : 
: 0?w 
= D sin 2(¢— 3) yt ; 
1 0?w 
M, = —+{D, + D.— (D,— Dj cos2g—- O } 
1 0? w (8) 
ae {D, + Ds — (Dz — Ds) cos2 (gy — #) Oy? . 
; 0?w 
+ D,sin2 (vy — #) 5) ay: 
1 072w 02w 
Mag = —>{ (Di — D2) Sox + (Do — Ds) Gar | sin2(p — 9) | 
0? w 
+ Dy, cos2(p — 9) =~ re | 
where 
gew {( 2%’ cosé 2 y’ sind \2 ‘cos? sin?}\ / #3 Ne \ 
Ox2 =2A( wr) b2 i 2( a2 b2 )( rT) a — i), 
Xm {/ 2x’ sind 2 y’ cos} \2 sin?d eoste \/ #3 yr 
oy? = 24] at eS ] + 2{ a? ~ Ba heer + ah 
02w 2’ cos# 2x’ sin® \ /2 x’ cos® 29’ sind’ : 
Ox Ov =2A|( be a a? \( a2 b2 . 
a8 ay oe orheraal eae 
Shen Sowa) — 1)(q: — a2)sin2 9}. 


As an example, we consider an oak plate cut parallel to the grain of wood. 
The axis of y is taken in the direction of the fibres, then it becomes!) (h being 
assumed to be unity), 


D,= 10-19, D,= 6-687, D, = 22-62, D, = 6-667, 
from which 


ky = 0-3428 — 0:2279 i, 3 = 0-3428 + 0-2279 i. 


1) H. Horie, Ing.-Arch. 6, 8 (1935). 
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The maximum deflection is at the center of the plate, and from equation 
(6), we obtain 
¢ ed A |) Ee Gee hie 
max L, cos*# + 2 L, cos? sin? + L, sintd ~ (9) 


Taking the ratio of the major and minor axes as 3:1, the amounts of the 
maximum deflection for three different cases of = 0, 2/4, 2/2 are calculated 
from (9). They are 


=0-1755, 0:2542 and 0:-3636 (unit 10-2 ¢ b4). 


70) 
Om ax 


When a plate is clamped at the periphery, the maximum moment will be 
there, and so it becomes essential to find the values of the moments at the 
boundary, from equation (8) it becomes 


a 


Wee 4A {Di + Ds + (D, — D,) cos2 (p — 9)} 


cos# cosB sin? sinf \2 
eg 
+ {D, + D3 + (D, — Ds) cos2 (y — d)} 


< sin? cosB eee i 2 D,sin2(y — 9) 


a 


> ( cos® cosp os sin? ne) (San eore _ _cos@ sinf ) 


a b a b 
M, =—4A |{D, + Dz — (D, — D,) cos2(y — 9)} 
7 (<osheese i sing sme | 
{ (10) 
Br Date a (D, — Ds) cos2 (p — 3)} | 
‘sin? cosp - cosé sinf \2 : 
x ad — ~} 4+ 2D,sin2(y — 9) 
‘cos? cosB sin? sinB \ (sind cosp — cosé smh) 
, (eabeee + anbnt aban sant 
‘ cos® cosh sin? sinp \2 
Mee = —4A {2 — D,) las Gree SS 5 ) 
+ (Dy — Dy) (SES cosh _ Sek SE) sin 2 (9 — 9) 
yD sin® sinB \ 
4 2D, cos2 (p — 9) (<ostcore {| 
sin? cosB cos® ome) 
x (—S Ale b | , | 
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where 
in2 2absin 2p 2 (a* + b?) cos 2B — a® + B® 
siné p = Gay 6? — (a? — 6%) cosap’ “°"“? ~ att 62 — (a? — 6%) cos 2p © 


The bending moments at the periphery are calculated from equation (10) 


ae Ny uk 


and are shown in figures 1, 2, and 3. For the convenience of comparison, the — 


corresponding values for the isotropy are also shown by dotted lines in the same 
figures (Poisson’s ratio is assumed to be 0-3). 

In the foregoing calculations, we treated only the case of a uniform load, 
but when the distribution of load is expressed by a polynomial of x’, y’, a 
similar procedure can be used by assuming A in (6) as a polynomial of x’, y’. 


2. Supported Edge 


The analysis for this case is more complicated. We transform the Cartesian 
coordinates into the following systems of curvilinear coordinates 


xtiy =cecosh(«+7 8), 
x+ik,y=c'e™ cosh (a’ +B’ +7,), (11) 
xtik,y=c"e'’ cosh(x”+iB” +7 8,), | 


where c’, c”, 0, #”, B,, and f, are arbitrary constants. In the transformation 
(11), from the conditions that the curve «’= a9, represents the periphery of 
the plate and that 6 = f’ at the boundary, we have 


c cosh & cos # = c’ cosh ay cos # cos 8, + c’ sinh a sin # sin B,, 

c sinh %& sin # = — c’ cosha cos ® sin B, + c’ sinha sin ® cos f,, 

ky ¢ sinh % cos ® = c’ sinh %& cos ® sin B, — c’ cosh a sin ® cos f;, “a 
— ky ¢ cosha sin 3 = c’ sinha cos ® cos B, + c’ cosh a sin & sin p,, 


from which 


(k? — 1) sinh 2 a sin2 # tan? p, | 
+ 4 {ki (cosh* x sin? — sinh? %,» cos?) + cosh? a» cos? d (13) 
— sinh® x, sin? 3} tan B, — (Ai — 1) sinh 2 % sin2 & = 0. 


Solving the equation (13), we have two solutions for tanf,, but the differ- 


ence of the corresponding values of f, is 2/2, and so it is sufficient for further 
calculations if we take only one value of the two. 
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Substituting the one value of f, into (12), we can find a», 3’, and c’. 

In the same manner, we can determine fy, «5, #”, and c” so as to satisfy 
the cognate conditions. If we determine the constants in this way, «=, — 
a’ =a, and «”=«% represent respectively the periphery of the plate and 


moreover 
all ed (14) 
at the boundary. 


Taking the deflection of the plate in the expression 


w= Rifi(x +7k, vy) + fo(x +7 Re y)] | 
+ Cy x4 4+ Cy x? y+ Cy x’ y24+C, x’ y’? (15) 
+0, "4 Cy ea C, x ae Cee | 
and 
fi(w++7k,y)= D3: (A, +4A)) cosh2 n (a’ + 7 B’), 
a (16) 
fi(*+¢k,y) = a, (B, +7 Bi) cosh2 n(x” +76"), | 
we have 
Cli cosh (2 7 + 2) a’ : | 
ee A 4 Mies jy Cos(2m + 2) B 
ccd | 1 1 Tamed hd ee é 
\(2n+1)2n 5 2n(2"— 1) J cosh 4% a, cose f 
cosh (2 2 — 2) «’ ; 
+ (2 m — 1) (2 » — 2) COS |S Fs 2) 6’ 
rag isan 4! sinh(Z2”-+ 2)a’  . 7 ; 
4 ae n | (2n+ 2) (anti) (@n+2) 8 
1 1 3 : 
Vita alpan tae (2 = 1) | sinh 2 ma’ sin2 n 6’ (17) 


Siniv(2 7 — 2 


) 
“I (2n — eee sin (2 — 2) 6’ 


ce"? Wo cosh (2 ” +- 2) o” : I 
+23 B, |e Ew ay 008 (2% + 2) B 


n=2Z 


1 1 , 
— he eahawe th 2m Bw iy) Cosh 2 na" cos2 m B" 


cosh (2 n — 2) a” 


* (2% — 1) (2 — 2) cos (2 n — 2) pr) — 
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c”2 ies) 


D sinh(2m + 2)o” . : 
4 we. [+3 n ae 2) (2 n + cay sin (2 n+ 2) B 


n=2 


A 


u Wines 
ES baw 1) 2n - a sinh 2” «" sin2n B 


1 is = PD aus ‘ 17 
ie s&s nh (22 2) sna (2 Phe 2) "| | ( ) 


(2m — 1) (25-2) 
+ Cy x4 Cy 8 4 Cy a8 9 4 Cy x’ y+ Cy 
+ C, x24 C, x’ y’ + Cy’? + Cy. 
To satisfy the differential equation (1), we have 
+ 3{C, (D, — D, — D4) — C, D3} cos* sin 8 
+ {3 C, (Dy + D, + D, + Dj) + (Dz + Dj) (6 C, — 7C, + 6 C,)} 
Gos7 7) sin? ? + 3.{C,(D,—D, + D+ C, D,\sin® 8 cos t 


(18) 


—— 


+ {3 0; D,+3C,D,+C, (Dz + D,)}sintd = 2. 


From the condition that w vanishes at the boundary, we find 


Bec Cosh O41 Os ¢ sinh* xs cosh*a) 4-3 Cxc* sinh*a, 


Ge cosheo, — C, c*sinh*a, + C, c* cosh*a, — C, c* sinh*a, 


; A B 
He eos 2a 45 6 * cosh 2 a4 = 0, 


+ 4 C, ¢? cosh*a,5 + 4.C, c* sinh?a5+ 8 C, = 0, 
C, c* sinh a, cosh? a + Cy c*sinh® a» cosh a +.2-C, c? sinh %& cosh a, 
t 


f B: : 
= se e’* sinh 2 a — > ¢”* sinh 2 co 0 


C, c* cosh*% — C3 c4 sinh? a» cosh?a% + C, c* sinh*a 


ae Gz = AS) c’2 cosh 4 a + (= — a3) c”* cosh 4 cH =0, 


Cy c* sinh & cosh® a» — Cy ct sinh*a%» cosh % 


aN . " 
s) c”?2 sinh 4 a = 0, 
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(43 Pel 


c’2 cosh 6 % 
120 er) ° cosh 6 a 


By i "2 ae 
ae (ane -- et By = ies) c’* cosh 6 a = 0, 


we At He — ; 
(36 lag i68 | Cyn 


ME ore: 
ae Gen ee $a ) c’* sinh 6 a = 0, 


oe 0166. > 80 8 6 thes a4 Oh ee 2 0 6.8 0 Ok ow 8 eee Oe (a6 (kOe we ere elim a eee wt Meee) 66a) 


Y Aye Ms ( Lath ae ee es Ani — c’3 cosh 2 ” NX 
2n—1 2n—1 Tr be Fone eals Sn 

US a (i yl 1 Brix | 6"? cosh 2 n a _0 
+ (ty (sar +a) at Oy 1J 8n ‘ 


{ Ans ( b> eget ae Ley ee Anus | 6? sinh2n a 
8 2 


1 1 ' Bai ) c"*sinh2nog 
Nin (5 nat! 2ae y But aa 4 i] 8n 0, 
and from the condition that M, vanishes at the boundary, we have the relation 
4 


{(D, + D,) + (D, — D,) cos2 (vy — #)} 


oa 


= {(Ds + Ds),+ (Dz — D3) cos2 (vp — & 32 ar (20) 


+ 2 Dysin2(p— 8) sy = 0, 


where 


Me) lo “a 9) a cosh 2 Xy COs 2 20 cos 2 B= + sinh 2 a sin 2 0 sin 2 f — cos20 
cosh 2 % — cos 2B 3 


sin 2 (p — 8) sinh 2 a cos 2 # sin 2 B — cosh 2 a sin 2 # cos 2 8 + sin 2 a 


cosh 2% — cos2 B ; 


I 


0?w Bay , " 
oP a py (A, cosh 2 n a + B, cosh2 n a) cos2” B 


nN=2 


— = 


oo 
(ie Uy / . : © 
-- s (A, sinh 2 2% + By sinh 2 ma) sin2 n B + 6 C, c? cosh? ay cos? B 
nee 


+ 3 Cyc? cosh* a sin # cos & + Cy c® (sinh? ay + sin? 9) 
+ 3 C,c* sinh* a, sin d cos # + 6 C; sinh? a, sin? # + 2 G cos? # + 
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+ 2C, sin 3 cosd + 2 Cg sin? d + {6 C, c? cosh? a, cos? o 
+ 3 C, c? cosh*a%» sin # cos + Cz c? (cosh?a, sin? — sinh? a, cos? B) 
— 3 Cyc? sinh? % sin 3 cos } — 6 C; c? sinh? a, sin? }} cos 2 B 


+ (3 C, cos? 3 + 4 C, sind cos d + 3 C, sin? 9).c? sinh x gosh a, sin 2 B, 


a = yo (A, kj cosh2 mag + B,, k3 cosh2 » «{) cos2nB 


+ (A), ki sinh 2 no + By, k3 sinh 2 a) sin2 n B 


+ 6 Cy c? cosh? a sin? } — 3 C, c? cosh? a, sin # cos # 

+ C..67 (sinh*«,-,-cos*?) — 3 C, c*sinh*a, sin? cos? 

+ 6C,;c*sinh?«, cos23 + 2:C, sin? ? — 2C, sind cost + 2 C, cos? d 
A {6 Cj c? cosh? a» sin? d — 3 C, c? cosh?» sin # cos b 

+ C3 c® (cosh? % cos? # — sinh? %» sin? #) 

+ 3 Cyc? sinh2a, sin } cos # — 6 C, c? sinh?%» cos? 3} cos 2 B 

+ (3 Cy sin? & — 4 Cy sind cos # + 3 C, cos? B) c? sinh a cosh a sin 2 B, 


a (A, k, sinh2 naj + B, ky sinh2 a) sin2n B 
Ox OV = s 


x 


— J) (Al, ky cosh2 2 a + Bi ky cosh 2 2 a) cos2 nm B 


—+- a [= 6 C, ¢* cosh®a, sin 20 => Cyc? sim2 2 


+6C;c? sinh? a, sin 2 & — 2 C,sin2®+2C,cos2d 


+ 2C,sin2 & +3 C, c* cosh?a) cos2 # + 3 Cyc? sinh?x, cos2 3 
+ {— 6 C, c? cosh? a sin2 # + 3 Cy c* cosh* a cos 2 0 

+ C3 c2 cosh2 ay sin2 d — 3 C,c* sinh® x cos2 0 

— 6 C,c® sinh?» sin2 #}cos2 B 

= {= 3 Cy c? sinh % cosh % sin 2 # + 4 C, c? sinh x cosh %» cos 2 


+ 3 C,c? sinh a cosh % sin 2 3} sin 2 ]. 
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Mathematical difficulties are encountered for the determination of the — 


constants C,, Cy,..., Cy, As, By, A, By, ... to satisfy the equations (18), (19), 
and (20), under the condition that the functions /,; and /, converge, so we shall 
find an approximate solution by a method as shown in the following example. 


We consider the same plate as in the foregoing example, and treat first the . 


case of = 2/4, since the procedure is quite similar in other cases. As the 
first approximation, assuming 


A, = B,= 4, = Bi= A,= B,= Aj= By= ++ =0, 


we find from equations (18), (19), and (20) 


C, = 0-0000824 g, Cy = 0-000218 q, C, = — 0:000218 q 6, 
C, = 0-0000242 g, C; = 0-002361 g, C, = — 0-013551 q b2, 
C, = 0-001013 9, C, = — 0-001985 q b2, Cy, = 0-012291 g bt. 


I 


At the boundary, w vanishes but M, does not vanish, and the boundary 
condition is not fully satisfied. The residual bending moment at the boundary is 


nf ee 0074 cos 4 B + 0-0046 sin 4B 9b? 
aaa 1-250 — cos2p ; 


As the second approximation, assuming 


A, B, ; Ay ; By : A, = By 2 Aj - By = «= 0, 
we have 
C, = 0-0000935 g, C, = 0-000219 gq, C, = — 0-000224 q b?, 
Cz = 0-0000249 g, C,; = 0-002332 q, C, = — 0-013406 q 02, 


C, = 0-001060 g, Cy = — 0-002088 q b2, Cy = 0-011194 ¢ d4, 
Ay = (0-0004465 — 0-0003158 7) g b2, By = (0-0004465 + 0-0003158 4) g 2, 
A,= (—0-0000297 +- 0:0000198 1) q b?, B; = (—0-0000297 — 0:0000198 ) q b2, 
and the residual moment is 


0:0052 cos 6 B — 0: <j 
M. = cos 6 B — 0-0040 sin 6 B qb? 
1-250 — cos 2B 3 


In the same manner, we have the residual bending moments for the second 
approximation, in the cases of # = 0 and 2/2, 


0:0037 cos 6 B 


M, = ~ “NiS80.— cos 2 B q6*, and M,, Sr 


0:0078 cos 6 B 
1-250 — cos2 8B 


q b?, 
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respectively. The bending moments at the center of the plate, where 


C= i= f=  f == 8” o =o > 


are 
IME = i iD D id 
o Pat 1 2 — (D; — D,) )cos2 3} —— 
1 
2 {Dz Ds — (Dz — Ds) cosa 8) 8 —D, sin2 a a 
: 0?w 
My = — > {D1 + Dz + (D; = D,) cos2 9} | 
1 (21) 
— 7 1D. + D, + (D, — Ds) cos2 9} 55 + Dy aie yee oe ee 
: 0>w 0°w) . 
Mi, = 2 {(D, D>) ava 7 (D, — Ds) a sin2 & 
— Dy cos 2.9 ae | 
in which 
2 
sae = 2C,cos?# + 2C,sin# cos? + 2C,sin?3 + A, + Bo, 
2 
a =2C,sin?7? — 2C,sind cost + 2C, cos? = A, kh? — B, Re, 
07w ; : : : 
De eo Gem Ca cosz vo, ©. sin 2h aA, k= Boks. 


Substituting the numerical values of the constants for the second approx- 
imation into (21), we have the bending moments for ? = 7/4, as 


M,, = 0-4241 ¢ 62, M,=0-2764 qb, M,, = — 0-0956 q b®. 


a 


In the same manner, we have the bending moments for ? = 0 and z/2, as 
M,, = 0-4572 ¢ b?, M,=0-1607¢ 07, M,,=0, (®=0), 
: é mG 
M,, = 0-4088 q 62, M, = 0-3689 9%, M,,=0, (9 es 5). 


From the foregoing calculations, it may be seen that the residual moment 
is a smaller quantity compared with the maximum moment, and it decreases 
rapidly according as we take more terms in the series in (17), and so if we take 
only a few terms of the series we can find an approximate solution with sufficient 


accuracy. 
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Next, we shall find the deflections on the axes of x’ and y’, using the 


solution for the second approximation. On the axis of x’, we have 


Mee Wetec ee ay at § 
w= (42 + St) x +(e, - 4 (Se +e) 


| 


Toye (22) 
1 P 1 72 
+ {Ce Fa ea oe B,)| ales ot OY te (A,¢? + Boo 7h 
and on the axis of y’, it becomes 
+ A kin Be RN a i te B,k\) , 
ae C18) ( i “3 on) yiaen Cs ake a os os 
(23) 
1 5 Oe oe Lae se 
f Ic, = (Ag hice B, ¥)| PO, Saag ee Boe ae | 
:10°gb" The deflections on the axes of x’ 
18 Riwnicer and y’, calculated from equations (22) 
-n---- Clamped edge and (23), are shown in figure 4. For the 
axis ¥ath convenience of comparison, the deflec- 
S a tions for a plate with a clamped edge 
hai ei Be by are also shown by dotted lines in the 
S gh eh same figure. 
3 
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Zusammenfassung 


Der Verfasser lést das Biegeproblem fiir eine diinne elliptische Platte aus 
orthotropem Material unter einer gleichmassigen Querlast. Er gibt Lésungen 
fiir zwei Hauptfalle an, namlich fiir eingespannten und aufgelegten Rand, und 
zwar unter der Voraussetzung, dass die Symmetrieachsen des Materials nicht mit 


den Ellipsenachsen zusammenfallen. 


Als Beispiel wird eine Eichenplatte numerisch durchgerechnet, die parallel 
zu den Fasern geschnitten ist und das Achsenverhialtnis 3:1 besitzt; die Ergeb- 


nisse werden mit denjenigen der isotropen Platte verglichen. 


(Received: October 17, 1951.) 
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Kurze Mitteilungen — Brief Reports — Communications bréves 


Die Knicklast des einseitis eingespannten, tangential gedriickten Stabes 


Von Max Beck, Ziirich!) 


1. Einleitung 


Damit, dass die Zahl der nichtkonservativen Stabilitatsprobleme viel kleiner 
ist als die der konservativen und dass sie in der Praxis relativ selten auftreten, 
lasst es sich wohl erklaren, dass sie meist ziemlich stiefmiitterlich behandelt werden. 
So ist zum Beispiel die Knicklast des einseitig eingespannten, tangential gedriick- 
ten Stabes (Figur 1) noch nie einwandfrei berechnet worden. Dieses Stabilitats- 
problem wird zwar mehrfach erwahnt?), wobei das Resultat, dass hier keine 
Knickgefahr bestehe, den Leser tiberrascht. Die vorliegende Arbeit wird nun 
aber zeigen, dass eine endliche Knicklast existiert. 

Die Anregung, dieses Stabilitatsproblem wieder aufzugreifen, gab eine Unter- 
suchung von H. ZIEGLER’), in der er zeigt, dass nichtkonservative Stabilitats- 
probleme nur mit dem kinetischen Stabilitatskriterium behandelt werden diirfen. 
Nach diesem ist die kritische die kleinste Last, unter der nicht mehr jede kleine 
Stérung zu einer Bewegung in der nachsten Umgebung der Gleichgewichtslage 
fiihrt. Das statische Kriterium, das nach der kleinsten Last fragt, ftir die eine 
nichttriviale Gleichgewichtslage existiert, kann hier zu ganz falschen Resultaten 
fiihren, wahrend es bei konservativen Problemen mit dem kinetischen Kriterium 
aquivalent ist. Diese Tatsache wurde bis jetzt nicht beachtet, und so ist auch der 
tangential gedriickte Stab nur statisch untersucht worden. 


2. Die Differentialgleichung 


Fiir das Studium der kleinen Bewegungen des urspriinglich vertikalen Stabes 
der Figur 2 diirfen die Funktion (4, ¢t), die die elastische Linie wiedergibt, und 
ihre Ableitungen als von erster Ordnung klein angenommen werden, und es ist 
erlaubt, von héherer Ordnung kleine Glieder zu vernachlassigen. Die Differential- 
gleichung der elastischen Linie lasst sich gleich wie in der Statik ableiten, wenn 

ei der Ermittlung der Beanspruchung der Stabquerschnitte auch die Tragheits- 
_ krafte beriicksichtigt werden: 


aT (8) = wl(E, t) oF dé. 
Die Dichte @ soll im folgenden langs des Stabes konstant sein, ebenso die Biege- 


steifigkeit «. Mit den Bezeichnungen von Figur 2 wird nach der Festigkeitslehre 
die Stabkriimmung y” bestimmt durch die Beziehung: 


1 
ay” = Py!) — ya] — HU 2) — f WE.) (E- *) oF ab. (1) 


1) Assistenz fiir technische Mechanik, ETH. es 

2) L. Corrarz, Eigenwertaufgaben mit technischen A nwendungen (Geest & Portig, Leipzig 1949), 
S. 41. — A. Pritcer, Stabilitdtsprobleme der Elastostatik (Springer, Berlin, Géttingen und Heidelberg 
1950), S. 217. 

3) H. Ziecier, Die Stabilitatskriterien der Elastomechantk, Ing.-Arch. 20, 49 (1952). 
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P’ ist in erster Naherung gleich der tangentialen Druckkraft P; in den Knick- 
fallen der Figuren 3a und 36 liefert auch die Lagerreaktion einen Beitrag zu H. 
Leitet man (1) zweimal nach x ab, so erhalt man fiir die kleinen Schwingungen 
des Stabes die Differentialgleichung 


en eer) (2 


Da sie H nicht enthalt, gilt sie auch fiir den durch eine axiale Kraft belasteten 
Stab. Die Lagerung der Stabe der Figuren 3a und 36 liefert die gleichen Rand- 
bedingungen fiir eine axiale wie fiir eine tangentiale Druckkraft, so dass die 
Eigenschwingungen dieselben sind. Die bekannten Eulerschen Knicklasten gelten 
hier also auch fiir die entsprechenden nichtkonservativen Knickfalle. 

Im freien Stabende verschwindet bei tangentialer Belastung sowohl die 
Querkraft Q als auch das Biegemoment, wahrend bei axialem Druck Q von Null 
verschieden ist. Die Lésungen von (2) miissen fiir den Stab der Figur 1 also die 
Randbedingungen 

(0) yi (O) ses (Eng ee 0 (3) 
erfiillen. 


3. Bestimmung der Knicklast; Diskussion 
In der partiellen Differentialgleichung (2) lassen sich mit dem Produktansatz 
y(%, t) = Y(*) TQ) 


die Variablen trennen, wobei man mit einem Exponentialansatz sowohl fiir Y 
wie T die allgemeine Lésung findet: 


y (x, t) = Ay eax eiat, 
Damit wird (2) zur charakteristischen Gleichung: 
adAt*+ PA?— oF wm’=0., 
Mit den Abkiirzungen 
; f 7) ee 3 aan ge RS i * ie 
Prfsigsie FA a Me \ ea: > a Aes veers PtP 


, 
a 2 


lautet die allgemeine Lésung fiir Y: 

Y(¥) = A cosh(A, *) + B sinh(A, ¥) + C cos(da, x) + Dsin(A, x) . 
Zur Ermittlung der vier Konstanten A, B, C und D geben die Randbedingungen 
(3) vier homogene lineare Gleichungen, die nur dann eine nichttriviale Lésung 


besitzen, wenn sie voneinander linear abhangig sind. Das ist der Fall, wenn die 
Eigenwertgleichung 


0= f(w%, P) 
= (2am? + p*) + 2a w* cosh(A, 1) cos(Ayl) + Va w® p sinh(A, 2) sin(A, 2) J (4) 
erfiillt ist. Sie ordnet jeder Belastung P unendlich viele diskrete Eigenfrequenzen zu. 
Nimmt P von Null an zu, so ist die kritische Belastung dann erreicht, wenn 


zum erstenmal w = 0 die Gleichung (4) erfiillt oder wenn (4) fiir eine Doppel- 
wurzel besitzt. Im ersten Fall existiert eine Lésung 


y(¥, t) = y(x), 


ee 


ee EN 


2 tee od 
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Pig. 1 
Der einseitig eingespannte, im freien Ende 
tangential gedriickte Stab. 


ee fo p)=0 
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Fig. 2 


Zur Berechnung des Biegemoments im 
Stabquerschnitt F. 


Fig. 3 
Zwei weitere Knickfalle 
unter tangentialem Druck. 


fog 20? 


0 225 


Fig. 4 


Bestimmung der Knicklast aus der Eigenwertgleichung. 


228 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


also eine nichttriviale Gleichgewichtslage. Die so gefundene Knicklast hatte 
auch statisch ermittelt werden kénnen. Im Falle einer Doppelwurzel kann y die 
Form 

y(x, t) = Y(*) tsin(o@ 2) 
annehmen. Die nichtkonservative Druckkraft leistet bei dieser Bewegung Arbeit, 
so dass die Amplitude der Schwingung immer grésser wird. Eine solche Knicklast 
kann nur mit dem kinetischen Kriterium gefunden werden. 

Gleichung (4) besitzt fiir kein endliches P eine Lésung mit w = 0. In einem 
Koordinatensystem mit m? als Abszisse und # als Ordinate (Figur 4) wird (4) 
durch eine unendliche Zahl von Kurvenasten wiedergegeben, die aus dem vierten 
Quadranten aufsteigen, im ersten Quadranten ein Maximum erreichen und wie- 
der in den vierten zuriickfallen. Das der Ordinatenachse nachste Maximum, das 
in Figur 4 wiedergegeben ist, besitzt den kleinsten p-Wert und bestimmt damit 
die Knicklast. Die numerische Auswertung von (4) ergibt: 

a 
— 20,05 Rr . 
Beim einseitig eingespannten, im freien Ende tangential gedriickten Stab ist also 
die Knickkraft achtmal grésser als die entsprechende Eulersche Knicklast. Lasst 
man die axiale Druckkraft tiber gemischte Belastungen stetig in eine tangentiale 
iibergehen, so wird dieser hohe Wert verstandlich. 

Das hier gefundene Resultat gilt vorlaufig nur fiir einen idealen, ddmpfungs- 
freien Stab. H. ZIEGLER?) hat gezeigt, dass bei nichtkonservativen Systemen die 
geringste Dampfung die Knicklast unter Umstanden erheblich verkleinern kann. 
Es wird beim vorliegenden Problem aber nicht ganz einfach sein, einen physi- 
kalisch einwandfreien Ansatz fiir die Dampfung zu finden, der die mathematischen 
Schwierigkeiten nicht betrachtlich erhGht. 


E, 


i 


Summary 


The author calculates the buckling load of a built-in free rod subjected to 
tangential end-thrust. As this stability problem is nonconservative, it cannot be 
solved by the usual static methods. Investigation of the free vibrations under P 
yields the buckling load P;, = 20,05 «/l? («: flexual rigidity, /: length). (Damping 
forces are not taken into account; they might diminish the actual critical load). 


(Eingegangen: 28. 1. 1952.) 


Bergmans Integraloperator erster Art und Riemannsche Funktion 
Von PETER HEnrici, Washington, D. C.?) 
Zur Integration linearer partieller Differentialgleichungen zweiter Ordnung 


vom elliptischen Typus hat St. BERGMAN den Begriff der vollstandigen L6- 
sungssysteme entwickelt*), Darunter werden Systeme von Partikularlésungen 


1) H. Ziecier, Die Stabilitdtskriterien der Elastomechanik, Ing.-Arch. 20, 49 (1952). 

2) American University. Die folgenden Untersuchungen wurden gemass einem Vertrage zwi- 
schen der American University und dem National Bureau of Standards, Washington, mit Férderung 
des Flight Research Laboratory, U.S. Air Force, ausgefiihrt. 


3) Vgl. St. BerGMAN, The Kernel Function and Conformal Mapping, Amer. Math. Soc. 1950, 
S. 112 ff., sowie die dort angegebene Literatur. 
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der vorgelegten Differentialgleichung verstanden mit der Eigenschaft, dass sich 
jede regulare Losung der Gleichung durch eine endliche Anzahl dieser Partikular- 
lésungen beliebig genau approximieren lasst. 

Die Existenz solcher vollstandiger Lésungssysteme kann durch die Betrach- 
tung gewisser Minimalprobleme bewiesen werden. Dieser Beweis enthalt jedoch 
keine explizite Vorschrift zur effektiven Konstruktion der Systeme. BERGMAN 
behebt diesen Nachteil im Falle eines einfach zusammenhangenden Grundge- 
bietes G, indem er gewisse lineare Integraloperatoren @Q[/] angibt, die die Ge- 
samtheit der in © regularen Funktionen einer komplexen Variablen auf die 
Gesamtheit der in © regularen reellen Lésungen der vorgelegten Differential- 
gleichung abbilden. Aus der Vollstandigkeit der Funktionen 2” (n = 0, 1, 2, ...) 
fiir die analytischen Funktionen einer Variablen folgt dann die Vollstandigkeit 
der Funktionen Q[2z”) fiir die Lésungen der Differentialgleichung. 

Insbesondere definiert BERGMAN den Integraloperator erster Art durch die 
folgende Eigenschaft: Das Urbild einer Lésung u(x, y) der Differentialgleichung 
sel, wenn 

BS seat Wy Be = 87 0 anal GAGA SCs 


gesetzt wird, gegeben durch U(z, 0). Fiir diesen Operator erster Art werden ver- 
schiedene Reihendarstellungen gegeben. 

Im folgenden wird darauf hingewiesen, dass dieser Operator auf das engste 
mit der Riemannschen Integrationsmethode hyperbolischer Differentialgleichun- 
gen zusammenhangt und sich ausserdem in einfacher Weise durch die Grundl6- 
sung der betrachteten Differentialgleichung ausdriicken lasst. 

Die gegebene Differentialgleichung sei 

02u 07u Ou Ou 


= + —— Bie 0. 1 
Tu] = yn toy? eee alee ae ) (1) 


Wir setzen mit BERGMAN voraus, dass a, b und c ganze analytische Funktionen 
des Punktes P = (x, y) sind. Nach klassischen Satzen existiert zu (1) fiir jeden 
Punkt QO = (&, 7) eine Grundlosung von der Form 


HPO) = 7270) log PQ 4 s(P, 0)  PQO=— Ve Shit aty— 027 (2) 
welche als Funktion von P der zu (1) adjungierten Differentialgleichung 
O20 s Gani 2U: O(a v) 0(b v) 


=- ae 0 3 
sO Oy? Ox Ov oe (3) 


mv] 


geniigt. Mit den Koeffizienten von (1) sind auch vy und s ganze analytische Funk- 
tionen von P und Q; ferner ist 

(P,P) =1. (4) 

Mit Hilfe der Grundlésung lasst sich jede in einem elementaren Gebiet der 

reellen (x, y)-Ebene zweimal stetig differenzierbare Losung von (1) darstellen 


durch i ‘ 
U()— oa i {fe a = ae + au 0| cos (”, &) | 

(5) 
| le i yi i ley o| cos (”, n)| dsg , | 


wo u = u(Q) und v = v(Q, P) gesetzt ist und der Punkt Q den Rand Xs) von G 
(mit der innern Normalen m) durchlauft. Aus (5) folgt, dass u(P) in &© analytisch 
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ist. Dariiber hinaus leistet (5) auch die analytische Fortsetzung von w in jenes 
Gebiet des komplexen (x, v)-Raumes, in welches die analytische Fortsetzung des 


logarithmischen Faktors log QP von v fiir alle Punkte Q auf 06 méglich ist. 
Fiihren wir die Variablen 
ge aetiy, C=E+in, 2 = e— ty, Cf ty (6) 


ein, so wird 
log PQ = log V(2— ¢) (@* = 0%), 


und wir erkennen, dass der Logarithmus nur fiir z€ 0@© oder z* € 06 singular 
wird. Damit gewinnen wir 

Satz I1): Die analytische Fortsetzung von u ist méglich nach allen Punkten 
(z, 2*) des kartesischen Produktes © X © des Gebietes © und seines honjugiert 
komplexen Gebietes ©). 

Es liegt nun nahe, die Transformation (6) schon in (1) vorzunehmen. Eine 
kurze Rechnung ergibt, dass, wenn 


(a, V) == Ui (as 28) al ee (2, ee eee 


und 
1 3 1 
gesetzt wird, (1) dann iibergeht in 
02?U 0U 0U ” 
LU] = 5, o8 | aoe See ee (7) 


D,E und F sind wieder ganze Funktionen von z und z*. Eine entsprechende 
Rechnung ergibt fiir (3) die zu (7) adjungierte Gleichung 
orFV 0(DV) O(EV) — 


TV Sede, a ee (8) 


(7) ist eine in z und z* hyperbolische Differentialgleichung. Beschranken wir 
uns fiir den Augenblick auf reelle Werte von z und z*, so kénnen wir sie mit Hilfe 
der Riemannschen Integrationsmethode integrieren. Insbesondere kénnen wir 
so das charakteristische Problem lésen, das heisst diejenige Lésung von (7) be- 
stimmen, die auf je einer Charakteristik der beiden Scharen [die in der Normal- 
form (7) durch z = const und z* = const gegeben sind] gegebene Werte annimmt. 
Verlegen wir den Ursprung des Koordinatensystems in den Schnittpunkt der 
Anfangscharakteristiken und schreiben wir 


U(z, 0) = f(z), U(0, 2*) = f*(z*)  [f(0) = f*(0)] (9) 


vor, so ist die Lésung’) des charakteristischen Problems gegeben durch 

1) Vel. St. BerGMAN, Zur Theorie der Funktionen, die eine lineare partielle Differentialgleichung 
bejriedigen, Recueil math. [nouy. s.] 2, 1169-1198 (1939). | 

*) G bedeutet jetzt das dem soeben betrachteten Gebiet © der reellen (x, y)-Ebene vermdge 
der Abbildung (*, vy) > z = x + 7 y entsprechende Gebiet der Gaussschen Ebene. 

3) R. Courant und D. Hirpert, Methoden der mathematischen Physik, Bd. 2 (Springer, Berlin 
1937), S. 316. 


ry 


——" 


Vol. III, 1952 Kurze Mitteilungen — Brief Reports — Communications bréves 231 


¢ 
O(, o*) = FOO, 0;5,6%) = 10) + f Fle, 0; 6%) [5h + Be, 0) f) a: 


au 


0 - 
+ / Fo z*: C, C*) io + D(0, 2*) f*| dz* 
0 


II 


Qf, f). 


Hier bedeutet F(z, z*; ¢, ¢*) die sogenannte Riemannsche Funktion der Glei- 
chung (7), die den folgenden Bedingungen geniigt: 


(a) M([F] = 0 in bezug auf z und z*; 


(b) ad (Zen omtere bn GN Et (zerCtenG cee) == Orin 2 und GF: 


Me 6) Osha) D(C, 2) P(E 2" 5 6, C*)==00 inet, ¢ und ¢*; 
Kel een Gem cn Garey) =—aleamee umd CF 


Diese Riemannsche Funktion ist aber, wie schon HADAMARD!) bemerkt hat, 
identisch mit dem gemass (6) transformierten Koeffizienten des Logarithmus in 
der Grundlésung (2) von (1), das heisst, es gilt 


Ud earn ace, ae ah (See tO iG) ey (A eg (12) 


Der Beweis dieses iiberraschenden Sachverhaltes fliesst sehr einfach aus dem 
Umstand, dass kraft (3) und (6) V(z, 2*; ¢, €*) = v(#, v; €, 7) der Differential- 
gleichung (8) geniigt: 


Bi MERIC 2 SEF Sa 
OR 


1 
et Fy , Pare E R| + regulare Funktion = 0. 


sa = DR 


Durch Vergleich der Singularitaten auf beiden Seiten und unter Beriicksichtigung 
von (4) folgen hieraus fiir R die Bedingungen (11). Wegen der (klassischen) ein- 
deutigen Bestimmtheit der Riemannschen Funktion folgt hieraus (12). Insbe- 
sondere ist damit auch gezeigt, dass die Riemannsche Funktion der Gleichung (7) 
eine ganze Funktion ist. Damit ergibt sich 


Satz II. Wahlt man f(z) und f*(2*) als in © bzw. 6 analytische Funktionen, so 
ist die durch (10) (mit F = R) definierte Funktion U(€, €*) in © x © analytisch. 


Speziell ist sie auch fiir ¢* = ¢, das heisst im reellen Gebiete ©, regular. 
Aus den Satzen I und II ergibt sich, dass die Formeln (9) und (10) eine ein- 
eindeutige Abbildung der Menge der in & regularen Lésungen von (1) auf die Menge 


der Paare (f, /*) von in © bzw. © regularen Funktionen einer Variablen [mit 
#(0) = f*(0)] liefern. Da gemass (10) Q[f, f*] stetig und linear von f und /* abhangt, 
folgt damit auch die Vollstandigkeit des Systems 


fz 0), OO, <2") Veet 0, 12a. 5 0 med, <2 


1) Courant-Hizs_ert, |. c., S. 431. 


232 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


Wenn die Grundlésung von (1) bekannt ist — und das ist bei vielen physikalischen 
Problemen der Fall —, so kann der Operator (10) in geschlossener Form dargestellt 
werden. Ist u(x, y) fiir reelle x und y reell, so ist U(0, z*) = U(z*, 0), und (9) wird 
mit der Umkehrformel fiir den Operator erster Art identisch. Daraus folgt die 
Identitat von (10) mit diesem Operator’). 

Es ist sehr instruktiv, die geschilderten Verhaltnisse am Beispiel der Wellen- 
gleichung zu verfolgen. Wir hoffen, diese Diskussion an anderer Stelle geben zu 
k6nnen. 


Summary 


ST. BERGMAN transforms by means of the so-called Integral Operator of the 
first kind analytic functions of one complex variable into solutions of a given 
partial differential equation of elliptic type. It is shown in the present paper 
that this transformation can be interpreted as the solution of a certain hyper- 
bolic equation by RIEMANN’s method. Furthermore, BERGMAN’s operator is ex- 
pressed in terms of the fundamental solution of the considered elliptic equation. 


(Eingegangen: 26. 2. 1952.) 


Zum Ergodenproblem 


Von Rupotr Kurtu, Bern?) 


Ubersicht 


Auf die Frage nach der Ersetzbarkeit des Zeitmittels einer Phasenfunktion 
durch ihr Phasenmittel (§ 1) wird eine elementare, aber fiir die konkreten An- 
wendungen noch am ehesten brauchbare Antwort gegeben (§§ 2 und 3).e Der 
Erwartungswert der mittleren Verweilzeit ist dem durchlaufenen Volumen zwar 
proportional, aber die individuellen Abweichungen kénnen relativ gross sein (§ 4). 


§ 1. Die Fragestellung 


Es liege ein mechanisches System von endlich vielen Freiheitsgraden vor. 
Seine Hamilton-Funktion sei von der Zeit ¢ unabhangig. Die Hamiltonschen 
Gleichungen mégen bei beliebig vorgegebenen Anfangswerten fiir alle ¢ eine ein- 
deutige, stetig differenzierbare Lésung besitzen. Sie definieren dann im Phasen- 
raume /' eine einparametrige Transformationsgruppe — anschaulich gesprochen: 
eine stationare Strémung. J schliesslich sei eine dieser Gruppe gegeniiber inva- 
riante Teilmenge des Phasenraumes; sie sei (LEBESGUE) messbar und von be- 
schranktem positivem Masse. 

1) BERGMANS Operator zweiter Art (vgl. etwa Sr. BercMan, The Approximation of Functions 
Satisfying a Linear Partial Differential Equation, Duke math. J. 6,537-561 [1940]) kann in folgender 
Weise als Modifikation des Operators erster Art gedeutet werden: Die Menge der Paare (f, /*) wird 
nicht direkt auf die Menge (L) der Lésungen von (1), sondern zunachst in nichttrivialer Weise auf 
sich selbst und erst dann auf (L) abgebildet, was gelegentlich (zum Beispiel im Falle der Differential- 


gleichung der kompressiblen Flissigkeiten) rechnerische Vereinfachungen mit sich bringt. 
2) Astronomisches Institut der Universitat Bern. 
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Ubrigens geniigt es fiir das folgende, wenn die Gruppe nur fiir fast alle An- 
fangspunkte x¥«€ J definiert ist. 

fo(¥), x € J , sei nun irgendeine in J definierte und summable Funktion (eine 
«Phasenfunktion»). Bedeutet x(t) das von der Gruppe zur Zeit ¢ erzeugte Bild 
des Anfangspunktes x und f(x, ¢) die Funktion fy [(¢)], so lautet eine Folgerung 
aus dem Liouvilleschen Satze: f(x, t) ist tiber J summierbar, und es ist 


| itz, t) dx =| f,(x) dx. 
J J 
Darin ist enthalten: Die Strémung ist inkompressibel. 
Ferner gilt der Satz von BIRKHOFF: Das «Zeitmittel» 


t 
pie vy ds 


existiert fiir fast alle ve J. 

Vielfach ist es tiblich, J auf einer Flache konstanter Energie zu wahlen. 
Begriffe und Satze sowie der folgende Satz tibertragen sich auf diesen Fall ohne 
weiteres bei Einfiihrung eines invarianten Flachenmasses. 

Die Ergodenhypothese besagt nun: Das Zeitmittel ist fiir fast alle ve J 


gleich dem « Phasenmittel » fe definiert durch 


Es ist bekannt: Die Hypothese gilt dann und nur dann, wenn J «metrisch un- 
zerlegbar» (metrisch transitiv) ist. Metrische Unzerlegbarkeit bedeutet: Fiir jede 
messbare invariante Teilmenge J’C J ist entweder J’ oder J — J’ Nullmenge. 

Indessen kennt man kein allgemeines Kriterium fiir metrische Unzerlegbar- 
keit. So scheint durch die neue Formulierung der zu beweisenden Behauptung 
wenigstens bis jetzt kaum etwas gewonnen zu sein. 


§ 2. Zur Antwort 


Zunachst: Die Ergodenhypothese gilt sicherlich nicht allgemeim — und vor 
allem nicht in den fiir die Anwendung wichtigen Fallen: Fiir die meisten «ver- 
niinftigen» mechanischen Systeme gelten die Satze von der Erhaltung des 
Impulses und des Drehimpulses, und fiir deren Komponenten trifft die Behaup- 
tung = f) im allgemeinen nicht zu. Die entsprechenden invarianten Teilmengen 
sind metrisch zerlegbar. 

Vergegenwartigen wir uns den Zweck der Ergodenhypothese: Es soll die 
Ersetzung des in der statistischen Mechanik benotigten, gewoOhnlich aber unbe- 
kannten Zeitmittels durch das leichter zu berechnende Phasenmittel gerechtfertigt 
werden. In der erstrebten Weise ist das allgemein also nicht méglich. Aber in 
anderer. Der folgende Satz gibt hiertiber Auskunft: 

Satz: 7. Der Erwartungswert des Zeitmittels (das ist sein Phasenmittel) ist gleich 
dem Phasenmittel der Phasenfunktion selbst. 


3. Diese Abschitzung kann allgemein nicht mehr verscharft werden. 
In Zeichen: 


ti, fa fe fiir alle f,; 


2. (f— fe < (fy — f,)? fiir alle fy; 


Sey road) 2 aay 
Der Satz erméglicht also die Anwendung der naturwissenschaftlichen Methode: 
Ersetzung eines unbekannten Individualwertes durch den bekannten Erwar- 
tungswert und Abschatzung des mittleren Fehlers. 
Beweis: Alle hingeschriebenen Mittelwerte existieren auf Grund der Voraus- 
setzungen und den Satzen von LiouvILLE und BirKHOFF zufolge. Es ist nun 


0 e: fy)? fiir mindestens ein f, . 


ie iP=f=h=fo 
nach dem Liouvilleschen Satze. 


(f- At = Fit =P te 


nach 1. Mittels der Schwarzschen Ungleichung folgt weiter 


is) 


Psp, 


pePp=R, 
(f—A* SR — fo? = (fo — fo)? - 
3. Man wahle fiir f, ein zeitunabhangiges erstes Integral der Hamiltonschen 
Gleichungen, etwa eine Komponente des Drehimpulses. 


§ 3. Korollar 


Der bewiesene Satz gilt wértlich auch fiir Zeitmittel iiber beliebige beschrankte 
Zeitintervalle. Ferner auch fiir Phasenmittel, die mit irgendeiner in J definierten 
positiven summierbaren invarianten Belegungsfunktion u(x) gebildet werden; ; 
fiir die Invarianz ist hinreichend und unerlasslich, dass uw allein von zeitunabhan- ; 
gigen ersten Integralen der Hamiltonschen Gleichungen abhangt. ' 

Die erste Bemerkung schrankt die Bedeutung des Birkhoffschen Satzes fiir — 
die Anwendungen ein. Die zweite kann wichtig werden, wenn es nicht gelingt, 
eine in Betracht kommende invariante Teilmenge J Cc J’ mit beschrinktem : 
positivem Masse wirklich zu bestimmen. Dann suche man mittels einer geeigneten _ 


Belegungsfunktion dem ganzen Phasenraume J’ ein endliches Mass zu geben und 
Setze., = J", 
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2. Die mittlere Abweichung der Zeitmittel gegen ihren Erwartungswert ist hoch- 
stens gleich der mittleren Abweichung der Phasenfunktion gegen thr Phasenmittel. 


§ 4. Anwendung 


M sei eine messbare Teilmenge von J und p(x) ihre charakteristische Funk- 
tion: 
fire + 2d. 


(2) = {5 
?M Lays fir x*eJ—M. 
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bo 
Oo 
On 


Dann ist fiir fo(¥) = py,(%) 


{2 Mass von M My 
a) Mass von he p 


die Wahrscheinlichkeit, dass irgendein Punkt x € J zu M gehort. Ferner ist 


+ Zeitdauer, ftir die te T und x(t) « MW 
Tae aS ae re 


die relative Verweilzeit des auf der Bahnkurve durch ¥ laufenden Phasenpunktes 
«(t) in der Menge M;|T | bedeutet hierbei die Lange des bei der Zeitmittelung 
zugrunde gelegten Zeitintervalls 7. Der Satz in § 2 liefert fiir die relative Verweil- 
zeit unmittelbar die beiden Aussagen - 


t=p, @—7'Sp(i-p). 


Die Behauptung «t = fiir fast alle ¥ ¢« J» trifft nur in Sonderfallen zu, zum 
Beispiel bei Systemen mit nur einem Freiheitsgrade und fest vorgegebener 
Energie. Im allgemeinen ist sie durch die soeben gemachten Aussagen zu ersetzen. 

Die Fehlerabschatzung lasst sich nicht ohne zusatzliche Voraussetzungen ver- 
scharfen, wie das Beispiel des harmonischen Oszillators zeigt: Fiir J und M 
wahlen wir das Innere je einer Bahnellipse, und zwar so, dass 


1 
Ml = JI. 
Dann ist 
t= yy (*) + Pf fiir alle ve J, 
BENS ae Sa 
T= De = * 
S 1 
[z= 2 = 2° 
a 1 Ree he ae 
Vice ete ahessnien 4 ae) eae, VT) 


LITERATURVERZEICHNIS 


A.J. Kuincuin, Mathematical Foundations of Statistical Mechanics (Verlag, 
New York 1949). 

E. Horr, Ergodentheorie (Verlag, Berlin 1937). Hier auch weitere Literatur- 
hinweise. 


Summary 


Let J be a measurable invariant part of the phase space of a given mechanical 
system, and let J habe the positive finite measure | J |. For f the following 
elementary statements are proved: (1) The mathematical expectation of the 
time average of any phase function coincides with its phase average. (2) The 
dispersion of the time average is not larger than that of the phase function. 
(3) The mathematical expectation of the time t of running through a measur- 


able set MC J equals p= |M|/|J|. (4) The dispersion of t is < //p (1 — ). 
(5) Generally the limits (2) and (4) cannot be reduced. 


(Eingegangen: 15. 1. 1952.) 
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Varia — Miscellaneous — Divers 


International Union of Pure and Applied Physics 


Under the patronage of the International Union of Pure and Applied Physics 
the two following meetings will be held during September 1952: 

‘““B and y Radioactivity”? at Amsterdam, from September 1st to 6th incl. 

Subjects 
Beta desintegrations. Allowed and forbitten transitions; /-capture. 
Gamma transitions and internal conversion. 
Systematics of the beta and gamma decay and the shellmodel. 
Angular and polarization correlation. 
Instruments for beta and gamma spectroscopy. 
Energy calibration of beta and gamma spectroscopes. 
Measurements of beta particles of small energy. 
Coincidence method for the analysis of complex spectra and for the 
determination of correlations. 

(9). Experiments of nuclear recoil and the detection of the neutrino. 

Detailed information can be obtained from the Secretary, Prof. Dr. J. DE 
Boer, Institute fiir Theoretical Physics, Roetersstraat 1A, Amsterdam, Nether- 
land. 

‘Structure and Properties of Solid Surfaces” at Chicago, IIll., from 
September 29th to October 1st incl. 
Subjects 

(1) Structure of crystal surfaces. 

) The thermodynamics of crystal surfaces; surface free energy, entropy and 
tension. 
3) Atomic mobility and reaction on surfaces. 
4) Crystal growth from vapor, nelt, and solution. 
5) Electrical properties of surfaces of metals, semi-conductors and insulators. 
6) 
7) 


NOU fWNH 
Rs NS er ee a 


oO 


SS a ee 


Surfaces of non-crystalline solids, particularly glasses. 

Cohesion and friction. 
8) Influence of surface on deformation of crystals. 

Detailed information can be obtained from the Secretary of the Organizing 
Committee: C.S. Smit, Institute for the Study of Metals, Chicago 37, Ill., U.S.A. 
A. PERRIER, Swiss National Committee of Physics. 


International Contest 


The Institute for the Unity of Science is offering a prize of $500 for the best 
essay on the theme Mathematical Logic as a Tool of Analysis: Its Uses and 
Achievements in the Sciences and Philosophy. Two additional prizes of $200 each 
will be given for the next best two essays. It is an International Contest and is 
open to everyone. Essays must not exceed 25,000 words. They may be written 
in English, French or German and must be submitted before January 1st, 1953. 
Further information can be obtained from the Institute for the Unity of Science, 
American Academy of Arts and Sciences, 28 Newbury Street, Boston 16, Massa- 
chusetts. PHILipp G, FRANK, Harvard University. 
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3. Osterreichischer Mathematikerkongress, Salzburg, 
9. bis 14. September 1952 


Die Osterreichische Mathematische Gesellschaft veranstaltet anlasslich des 
3. Osterreichischen Mathematikerkongresses wiederum ein Internationales Mathe- 
matikertreffen in Salzburg, in der Zeit vom 9. bis 14. September 1952. 

Die wissenschaftliche Arbeit des Kongresses wird sich in die folgenden Sek- 
tionen gliedern: Analysis; Geometrie und Topologie; Algebra und Zahlentheorie; 
Angewandte Mathematik; Geschichte und Philosophie. Bereits haben sich iiber 
400 Vertreter aus folgenden Landern angemeldet: Algerien, Belgien, Brasilien, 
Danemark, Deutschland, England, Finnland, Frankreich, Griechenland, Italien, 
Jugoslawien, Kanada, Niederlande, Nigeria, Norwegen, Polen, Schweden, Spanien, 
Schweiz, Sitidafrika, Tschechoslowakei und den USA. 

Anmeldungen und Anfragen sind zu richten an das Sekretariat der Osterrei- 
chischen Mathematischen Gesellschaft, Technische Hochschule, Karlsplatz 13, 
Wien IV. R. INZINGER 


Schweizer Archiv ftir angewandte Wissenschaft und Technik 
Zeitschrift fiir das Gesamtgebiet der Werkstoffkunde und der Werkstoffpriifung 


Seit Beginn dieses Jahres ist das Schweizer Archiv in neuer Aufmachung er- 
schienen. Im Geleitwort zur ersten Nummer des neuen Jahrganges wird darauf 
hingewiesen, dass die Zeitschrift sich in Zukunft auf das Gebiet der Werkstoff- 
forschung und Werkstoffpriifung beschranken und konzentrieren will. Da dieses 
Gebiet der angewandten Wissenschaften sich in den letzten Jahren auch in der 
Schweiz stark entwickelt hat, ist es sehr zu begriissen, dass dieser Forschungs- 
zweig sein eigenes Publikationsmittel gefunden hat. Zu unserer lebhaften Genug- — 
tuung ist es damit automatisch auch zu einer besseren Abgrenzung der Sachge- 
biete mit der ZAMP gekommen, die sicherlich nur zum Nutzen beider Zeitschriften 
gereichen wird. 

Die Schriftleitung des Schweizer Archivs verbleibt auch weiterhin in den 
Handen von PD. Dr. H. STAGER; ihm zur Seite steht eine Redaktionskommission 
mit einem Redaktionsausschuss, dessen Prasidium Direktor A. IMHOF inne hat. 

Die Redaktion. 


Errata 


Sur la thermodynamique des processus irréversibles. Par KyriLLr 
Poporr (ZAMP 3, Fasc. 1, 42 [1952)). 


Au lieu de X,, des formules (4), (5), (10) lire — X;,, c’est-a-dire 
AS) = — > X, Orn, Ao — — Sees 
k k 


au lieu de j 
O(AS) = DX; 0%, , ANS) = 2 A Bp 
k R 
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et 
#° + #3" pe rae Et Ss) =) Ge (11) 
2 1 
au lieu de : 
sect Be tt Bike Se 67, Kg? Sy (11) 
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Vorlesungen iiber theoretische Physik, Bd. III: Elektrodynamik. Von 
A. SOMMERFELD (Dieterichsche Verlagsbuchhandlung, Wiesbaden 1948). 367 S., 
48 Abb.; DM. 15.-. 


Inhalt: I. Grundlagen und Grundbegriffe der Maxwellschen Elektrodynamik; 
II. Ableitung der Erscheinungen aus den Maxwellschen Gleichungen; III. Rela- 
tivitatstheorie und Elektronentheorie; IV. Die Maxwellsche Theorie fiir bewegte 
Kkd6rper und andere Erganzungen. 

Im I. Teil werden die Maxwellschen Gleichungen (in Integral- und Differential- 


form) axiomatisch an die Spitze gestellt und vor allem die prinzipiellen und ~ 


praktischen Vorziige des Giorgischen MaBsystems herausgearbeitet. 

Im II. Teil werden die klassischen Anwendungen der Maxwellschen Gleichun- 
gen von den einfachsten elektrostatischen Problemen bis zur Dipolstrahlung und 
den Drahtwellen durchgerechnet. 

Teil III bringt vor allem eine sehr klare Ableitung der speziellen Relativitats- 
theorie aus der Forderung der relativistischen Invarianz der Maxwellschen Glei- 
chungen, Teil IV neben der Elektrodynamik bewegter K6rper eine fast elemen- 
tare Behandlung wichtiger Satze der allgemeinen Relativitatstheorie. 

Den Schluss bilden 25 sehr instruktive Ubungsaufgaben aus allen vier Teilen 
des Buches. 

Weder der Name SOMMERFELD noch die beriihmte Buchserie seiner Vovrle- 
sungen tibey theoretische Physik bediirfen einer Empfehlung; beide sind schlechthin 
«klassisch». Der Praktiker wird vor allem den klaren, deduktiven Aufbau des 
Hauptteils: «Ableitung der Erscheinungen aus den Maxwellschen Gleichungen» 
sowie die konsequente Anwendung des Giorgischen Mafsystems schatzen; der 
Theoretiker dagegen wird hauptsachlich die glanzende Ableitung der speziellen 
Relativitatstheorie aus der Invarianz der Maxwellschen Gleichungen mit Begei- 
sterung lesen. Das Buch ist weder textlich noch drucktechnisch als Nachschlage- 
werk geeignet; es ist ein klassisches Lehrbuch, sicher das beste, das je iiber dieses 
Gebiet geschrieben worden ist. Es stellt an den Leser hohe Anforderungen, 
bringt ihm aber grossen Gewinn. P. Scherver. 


Randwertprobleme und andere Anwendungsgebiete der héheren 
Analysis fiir Physiker, Mathematiker und Ingenieure. Von F. ScHwank 
(B.G. Teubner, Leipzig 1951), 406 S., 147 Abb.; DM. 22.50. 


Das Buch kann den Ingenieuren und Physikern, die in méglichst einfacher 
Weise in die Behandlung von Randwertproblemen und von verwandten Auf- 
gaben eingefiihrt zu werden wiinschen, warm empfohlen werden. Der Verfasser 
hat sich bemuht, die mathematischen Ableitungen einfach und ausfiihrlich dar- 
zustellen; auf Existenz- und Konvergenzbeweise hat er im allgemeinen verzichtet. 


on Ohl tanta 
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Nach der mathematischen Seite hin setzt die Lektiire des Buches nur die Kenntnis 
der Grundlagen der Differential- und Integralrechnung voraus; bei der Behand- 
lung von physikalischen und technischen Beispielen wird die Kenntnis der betref- 
fenden Grundgesetze ebenfalls vorausgesetzt. 

Das erste Kapitel gibt eine Einfiihrung in die Darstellung und Liésung der 
Randwertprobleme der schwingenden Saite und des schwingenden Stabes. Im 
zweiten Kapitel ist die Theorie der analytischen Funktionen einer komplexen 
Variablen bis zu den Anwendungen des Cauchyschen Integralsatzes entwickelt. 
Ein drittes Kapitel behandelt an Hand konkreter Aufgaben aus der Potential- 
theorie, der Theorie der Membran- und Plattenschwingungen und der Theorie 
der Warmeleitung die wichtigsten vorkommenden Differentialgleichungen und 
ihre Losungen (Legendresche Polynome, Kugelfunktionen, Besselsche Funktionen 
usw.). Das vierte Kapitel fiihrt in die Theorie der Integralgleichungen und ihrer 
Anwendungen ein, das fiinfte in die Variationsrechnung und das sechste in die 
Theorie und Praxis der Differenzengleichungen. Literaturangaben sind am Ende 
jedes Kapitels gegeben; in einem Anhang findet der Leser eine kurze Rekapitu- 
lation der Hauptsatze tiber Determinanten, lineare Gleichungssysteme, Kurven- 
integrale, unendliche Reihen, lneare Differentialgleichungen, die im Buch zur 
Verwendung kommen. M. Plancherel 


Nonlinear Vibrations in Mechanical and Electrical Systems. By 
J. J. STOKER (Interscience Publishers, Inc., New York 1950), 273 pp., $5.—. 


Das Buch ist aus Vorlesungen und Seminarien des Verfassers an der New York 
University entstanden und stellt den zweiten Band der von ihm zusammen mit 
H. Bour und R. Courant herausgegebenen Reihe «Pure and Applied Mathe- 
matics» dar. Es ist als Lehrbuch fiir Ingenieure und Physiker konzipiert und fiihrt 
in vorbildlich klarer Weise sowie in standiger Fiihlung mit praktischen Beispielen 
in die Theorie der nichtlinearen Schwingungen mit einem Freiheitsgrad ein. Was 
das Werk vor anderen Publikationen auf diesem Gebiet auszeichnet, ist die Kunst 
des Verfassers, bei aller Fiille des Stoffes die Akzente auf das Wesentliche zu 
setzen, ferner die strenge Systematik, mit der die Differentialgleichung 


mk + p(#) + f(x) =F cosat 


nach ihren praktischen Erscheinungsformen zergliedert und im Hinblick auf die 
esonderen Merkmale ihrer Lésungen diskutiert wird. 

Die einzelnen Kapitel umfassen: I. lineare Schwingungen (p=c%x,f=k 4), 
II. freie, ungedampfte nichtlineare Schwingungen (y = 0, F = 0), III. freie, ge- 
dampfte Schwingungen (F = 0), erzwungene Schwingungen mit IV. nichtlinearer 
Riickstellkraft (y= c #) und V. nichtlinearer Dampfung (f= # %) sowie VI. 
quasilineare Schwingungen [py = p(t): %, f= q(t): ¥, F = 0); eine Reihe von Exi- 
stenz- und Eindeutigkeitsbeweisen sind im Anhang zusammengefaBt. H. Ziegler 


Handbuch der Laplace-Transformation, Bd. I, Theorie der Laplace- 
Transformation. Von G. DortscH (Verlag Birkhauser, Basel 1950)\poSilssos, 
40 Abb.; sFr. 78.-. 

Es handelt sich um die zu einem Handbuch erweiterte und verbesserte Aus- 
gabe des 1937 erschienenen Werkes Theorie und Anwendungen dev Laplace- Tvans- 
formation des gleichen Verfassers. Der vorliegende erste Band korrespondiert 
dabei im wesentlichen mit den Teilen I, II, und III der alten Ausgabe, beschrankt 
sich aber auf die Theorie, die indessen erschépfend behandelt wird. Insbesondere 
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sind neben anderem die Parsevalsche Gleichung, die Erweiterung der Laplace- 
Transformation durch arithmetische Mittelbildung und die Umkehrung der 
Transformation durch Reihenentwicklung in der neuen Ausgabe weit ausfuhr- 
licher behandelt. Die Anwendungen werden gesamthaft im zweiten Band zu 
finden sein. 

Inhaltsiibersicht von Band 1: 1. Teil: Grundlegende analytische und funk- 
tionentheoretische Eigenschaften der Laplace-Transformation. Hier werden unter 
anderem behandelt: Konvergenzeigenschaften, Abbildung der Differentiation 
und Integration, Faltung; ferner auch Verhalten der Bildfunktion bei Annahe- 
rung an s = oo. 2. Teil: Umkehrung der Fourier- und Laplace-Transformation, die 
Parsevalsche Gleichung, Umkehrung durch Reihenentwicklung, Darstellungs- 
problem. 3. Teil: Die Cesaroschen arithmetischen Mittel des Laplace-Integrals 
und die 2X-Transformation. 4. Teil: Die Laplace-Transformation spezieller Klas- 


nal gmt 


sen von Funktionen (insbesondere von ganzen Funktionen vom Exponential-_ 


typus). 5. Teil: Abelsche und Taubersche Satze (Theoretische Grundlagen der 
Asymptotik). Am Schlu8 folgen ein Anhang mit Hilfssatzen und historischen 
Bemerkungen sowie ein ausfiihrliches Literaturverzeichnis. 

Dieses umfassende Werk wird der Theorie und Praxis gleichermaBen wert- 


volle Dienste leisten und kann als Nachschlagewerk fiir alle theoretischen Fragen — 


liber die Laplace-Transformation warmstens empfohlen werden. 


H. Rutishauser 


An Essay toward a Unified Theory of Special Functions, Based upon — 


the Functional Equation 
6) 
Oz F(z, a) = F(z,a@ Siz 1). 


By C.TRUESDELL (Princeton University Press, 1948). 188 pp.; $3.00. 


Les fonctions classiques de la physique mathématique (fonctions Beta et 


“ee 


Gamma, fonctions hypergéométriques, fonctions de Legendre, de Jacobi, de 


Gegenbauer, de Laguerre, d’ Hermite, de Weber, fonctions de Bessel et de Hankel, 


etc.) satisfont a un grand nombre de relations, disparates au premier abord. Ces © 
fonctions peuvent étre obtenues a partir de fonctions génératrices relativement 


simples ; elles s’expriment par suite par des intégrales de contour ou des intégrales 
définies et elles satisfont & des équations aux différences. 

L’auteur montre que toutes ces relations découlent du fait que ces fonctions 
sont des solutions d’équations aux différences du type général 


) 
jy 1 ®) = Aly, @) Hy, @) + Bly, a) fy, @ + 1) 


ou A et B sont des fonctions données et B # 0. Par une transformation con- 
venable, cette équation aux différences se raméne au type 


L’ouvrage est consacré a l’étude systématique de cette derniére équation. Cette 
étude permet non seulement de retrouver les formules connues en unifiant les 
théories des diverses fonctions citées plus haut; elle donne aussi plusieurs relations 
nouvelles et pose des problémes qui attendent encore leur solution. 


M. Plancherel 
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Geratehauanstalt Balzers 


Hochvakuum- 
Bedampfungsanlagen 


zur Herstellung dinnster Schichten im 
Vakuum. Besonders geeignet zur Optik- 
vergitung, Herstellung hochwertiger 
Spiegel und Beschattung von Prapara- 
ten fiir Elektronenmikroskopie. Dreistu- 
fige Ol-Diffusionspumpe hoher Saug- 
leistung, zweistufige Vorvakuumpumpe, 
VakuummeBgerate fiir 107° bis 1 mmHg. 
Unter Vakuum rotier. Drehvorrichtung 
zur Erzielung gleichmaBiger Schichten. 


Verkaufsorganisation 


Triih, Tauber & Co. AG., Ziirich 
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